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1. Introduction 

Polylogarithms and their multivariable generalizations [1, 2] play an equally important role 
in modern mathematics and in physics. In mathematics they occur for instance in connec- 
tion with algebraic i^-theory and mixed Tate motives, e.g. [3, 4, 5, 6, 7, 8], with Hilbert's 
third problem (on scissors congruences), e.g. [9, 10, 11], as volume functions for hyperbolic 
spaces, e.g. [12, 13, 18, 14, 15, 11, 16], and are also related to characteristic classes, e.g. [17], 
special values of L-functions in algebraic number theory, e.g. [18, 5, 7], algebraic cycles, 
e.g. [19, 20, 21] or, in the form of iterated integrals, in algebraic topology, e.g. [22, 23, 24]. 
In physics, the computation of higher order corrections to physical observables requires 
the analytical evaluation of Feynman integrals that can generally be expressed in terms of 
(special classes of) multiple polylogarithms, e.g. [25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 
36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57]. While 
in all of these applications it would be desirable to have a minimal spanning set — "basis 
functions" in physics parlance — for the polylogarithmic expressions involved in a given 
problem, it is well known that these latter functions satisfy various intricate functional 
equations among themselves, making the question of how to find a minimal spanning set 
very hard to answer in general. As a consequence, seemingly complicated results, say for a 
Feynman integral, may admit a much shorter analytic representation, the simplicity of the 
answer being hidden due to the existence of an abundance of functional equations among 
these functions. There is thus a strong interest for a better understanding of the functional 
equations among multiple polylogarithms, both from a formal mathematical standpoint 
and in view of practical applications in physics. 

A way to approach functional equations among (multiple) polylogarithms is provided 
by the so-called symbol map, a linear map that associates to each multiple polylogarithm 
of weight n an element in the n-fold tensor power of some vector space of one-forms. The 
virtue of the symbol map is that it captures to a good extent the main combinatorial and 
analytical properties of certain transcendental functions, and in particular it is expected 
that all functional equations among multiple polylogarithms are in the kernel of the symbol 
map. Loosely speaking, this means that a necessary condition for two expressions written in 
terms of multiple polylogarithms to be equal modulo functional equations is that they have 
the same symbol, a condition that is usually much easier to check than proving equality at 
the level of the functions. The inverse problem (sometimes called integration of a symbol) 
of finding a function whose symbol matches a given tensor satisfying a certain integrability 
condition is much harder and we know of no general algorithm to construct such a function. 

While special cases of the symbol map have been profitably used by mathematicians for 
over two decades (for example in connection with functional equations see e.g. refs. [5, 7, 
58, 59]), it has only very recently been introduced into physics in the context of the N = 4 
Super Yang-Mills (SYM) theory in ref. [60], where it was applied to greatly simplify the 
analytic expression for the two-loop six-point remainder function obtained in ref. [61, 62]. In 
the wake of that work, the symbol map has seen various applications, mostly in the context 
of N = 4 SYM. In particular, by now the symbols of all two-loop remainder functions are 
completely known [63] , while at three loops the symbols of the remainder functions for the 
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hexagon in general kinematics [64] and for the octagon in special kinematics [65] are known 
up to some free parameters that could not be fixed from general considerations. However, 
only in the latter octagon case an integrated form of the symbol is also known. Other 
approaches, aiming at the determination of the symbol of loop amplitudes by exploiting 
the operator product expansion in the collinear limit [66, 67] or the relationship between 
Feynman integrals and the volumes of polyhedra in non-euclidean spaces [68, 69], have 
also been considered. Furthermore, the symbol map was recently used to obtain compact 
analytic expressions for certain one-loop hexagon integrals in D = 6 dimensions [70, 71, 
72, 73]. More phenomenological applications, as for example in ref. [74], have also been 
considered. 

The aim of this paper is twofold: While the symbol map has already been extensively 
used in the N = 4 SYM community in physics, it seems still rather little known in other 
areas of physics in which the computation of Feynman integrals plays an important role. 
On the one hand, we therefore present a concise review on this topic, putting special 
emphasis on how to apply the symbol map to obtain simpler or shorter analytic results 
for functions arising from certain Feynman integrals. On the other hand, we believe that 
our work goes beyond the existing literature on the subject in various aspects. While so 
far the symbol of a transcendental function was defined recursively by considering iterated 
differentials, we introduce a simple diagrammatic rule that allows to directly read off the 
symbol from all possible "triangulations" of a certain decorated polygon associated to a 
multiple polylogarithm [21]. Furthermore, we also address the problem of how to integrate a 
symbol to a function by presenting an effective approach to construct a candidate spanning 
set of functions in terms of which the symbol might be integrated. 

The structure of the paper is as follows: In section 2 we give a short review of multiple 
polylogarithms and of their properties. In section 3 we review the main properties of the 
symbol map and we show how to obtain the symbol of a multiple polylogarithm as the 
weighted sum of all possible maximal dissections of a certain decorated polygon associated 
to the polylogarithm. In section 4 we give a short example of how to integrate a symbol 
following the approach introduced in refs. [60, 75], before generalizing this procedure to 
higher weights in section 5. In order to highlight remaining difficulties and ambiguities 
when trying to integrate to a function, we also give a family of non-trivial elements in the 
kernel of the symbol map. We illustrate these concepts in section 6 where we apply them 
to derive a spanning set up to weight four for a special class of multiple polylogarithms, 
the so-called harmonic polylogarithms [25]. The appendices contain a summary of the 
mathematical notions used throughout the paper, as well as some technical details and 
proofs left out in the main text. We also include an appendix with a collection of symbols 
for multiple polylogarithms up to weight four. 

Remark: The authors wish the reader to be aware that this paper contains the work 
of both physicists and mathematicians. As a consequence, it should be noted that the 
paper has been written to try to accommodate the language of both communities. We 
have tried to find a compromise in the level of details of the paper, and so while some 
arguments may go deeper than felt necessary by some readers, the text may be too sketchy 
at times for others. 
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2. Short review of multiple poly logarithms 



Definition. Multiple polylogarithms can be denned recursively, for n > 0, via the iterated 
integral [1, 2] 

f x dt 

G(a±, . . . ,a n ;x) = / G(a 2 , . . . , a n ; t) , (2.1) 

Jo t — ai 

with G(x) = G(;x) = 1, an exception being when x = in which case we put G(0) = 
(clearly any expression J ° . . . should be zero), and with at G C are chosen constants and x 
is a complex variable. In the following, we will also consider G(a±, . . . , a n ; x) to be functions 
of a±, . . . ,a n . In the special case where all the a,'s are zero, we define, using the obvious 
vector notation a n = (a, . . . , a), a G C, 



G(0 n ;x) = - log n x, (2.2) 



consistent with the case n = above. Note that, while in the Mathematics literature these 
functions appear already in the early 20th century in the works of Poincare and of Lappo- 
Danilevsky [76] as "hyperlogarithms" , as well as in the 1960's in Chen's work on iterated 
integrals (e.g., [22] ) 1 , in the physics literature these functions are often called Goncharov 
polylogarithms, due to the wealth of structure that the latter has established for them 
over the last 20 years. Throughout this paper, we follow the physics convention for the 
definition of the iterated integrals, which differs slightly from the mathematical one; e.g., in 
ref. [2], the function corresponding to G(a±, . . . , a n ; x) would be denoted 1(0; a n , . . . ,a±; x), 
i.e., with the reverse order of the aj but keeping the same variable x. 

We will refer to the vector a = (ai, . . . ,a n ) as the vector of singularities attached to 
the multiple polylogarithm and the number of elements n, counted with multiplicities, in 
that vector is called the weight of the multiple polylogarithm. 

Properties. We collect here a number of useful and well-known properties (cf. e.g. 
ref. [2, 8]). Iterated integrals form a shuffle algebra [78] (see appendix A for a short review of 
shuffle algebras) , which allows one to express the product of two multiple polylogarithms of 
weight ni and ri2 as a linear combination with integer coefficients of multiple polylogarithms 
of weight ni + ri2, via 

G(a±, . . . , a ni ; x) G(a ni +i, . . . , a ni + n2 ; x) = G^a^i), . . . , a .( rai+n2 ); x), (2 3) 

(reS(rii,ri2) 

where £(^1,712) denotes the set of all shuffles of n\ + n 2 elements, i.e., the subset of the 
symmetric group S ni + ri2 defined by (cf. ref. [22], eq. (1.5.6)) 

S(ni,n 2 ) = {a G S ni+n2 \ cr~ 1 (l) < . . . < cr~ 1 (ni) and cr~ 1 (ni + l) < . . . < cr~ 1 (ni+n 2 )} . 

(2.4) 

The algebraic properties of multiple polylogarithms imply that not all the G(a; x) for fixed 
x are independent, but that there are (polynomial) relations among them. In particular, we 



1 In a sense, they already made an appearance in Rummer's pioneering work [77] in 1840. 
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can reduce them, modulo products of lower weight functions, to functions whose rightmost 
index of all the vectors of singularities is non-zero (apart from objects of the form G(0 n ; x)), 
e.g., 

G(a, 0, 0; x) = G(0, 0; x) G(a; x) - G(0, 0, a; x) - G(0, a, 0; x) 

= G(0, 0; x) G(a; x) - G(0, 0, a; x) - (G(0, a; x) G(0; x) - 2G(0, 0, a; x)) (2.5) 
= G(0, 0; x) G(a; x) + G(0, 0, a; x) - G(0, a; x) G(0; x) , 

where the middle summand is of the desired form (and the remaining summands are prod- 
ucts). 

If the (rightmost) index a n of a is non-zero, then the function G(a; x) is invariant 
under a rescaling of all its arguments, i.e., for any k G C* we have 

G(ka; kx) = G(a;x) (a n ± 0) . (2.6) 

Furthermore, multiple polylogarithms satisfy the Holder convolution [79], i.e., whenever 
ai 1 and a n / 0, we have, Mp G C*, 

G(ai,...,a n ;l) = ^(-l) fc G^l-a fc ,...,l-ai;l-^ G ^a fc+ i, . . . , a n ; ^ . (2.7) 

Below in section 5 we will be particularly interested in the limiting case p — > oo of this 
identity, 

G(oi, . . . , a„; 1) = (-1)" G (1 - a n , . . . , 1 - a V , 1) . (2.8) 

Whenever they converge, multiple polylogarithms can equally well be represented [1] 
as multiple nested sums (e.g., for |xj| < 1) 

~rn nv n fe oo n k n k -l n 2 -l ni 

t ■ / A v ^ ^2 k \ ^ k \ ^ v ^ x-^ 

^mi,...,m k \ x li • • • , £fc) = 2^ „mi m 2 = 2_> ' ' ' ' 

I ley * * * /(/j, Ith, /ti 

ni<n 2 <---<n fc 1 z « rifc=l K n fc _x=l m=l 1 

(2.9) 

Note that we are using Goncharov's original summation convention [1]; other authors define 
Li mi ,..., mfc (xi, . . . ,Xk) using the reverse summation convention instead, i.e. ri\ > ■ ■ ■ > n k . 
The G and Li functions define in fact the same class of functions and are related by 



Li mi) ..., mfc (xi, . . . , Xfc) — (— 1) G mfcj ... jmi I , | , (2-10) 

(note the reverse order of the indices for G) where we used the notation 

G mi ,...,m k (h,...,t k ) = G((V^0, t u ... , (V^O, t k ; 1) . (2.11) 

mi— 1 vn k — \ 

It is possible to find closed expressions for (very few) special classes of multiple polylog- 
arithms, for arbitrary weight, in terms of classical polylogarithm functions, e.g., for a / 
we have 

G(0 n ; x) = - } log™ x, G(a n ; x) = - log™ (l - -) , 

n! n! V a) ^.12) 

-* / X \ -> / X \ 

G(0 n _i,a;x) = -Li n (-J , G(0„,a p ;x) = (-l) p S niP (-J , 

\ (X / \ CL / 
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where S np denotes the Nielsen poly logarithm [80]. Moreover, up to weight three, multiple 
polylogarithms are well-known to be expressible in terms of ordinary logarithms, dilog- 
arithms and trilogarithms (cf. ref. [81], §8.4.3, implicitly, as well as refs. [11, 15]). In 
particular, if a and b are non-zero and different, we find 

GCo.bs x) - IJ. - I*. (^) + log (l - f ) log (f^i) . (2.13) 

Aim. The aim of this paper is to present an algorithmic approach how to deal with — in 
fact rather to circumvent — the complicated functional equations that relate multiple poly- 
logarithms, and how to find, given a choice of certain singularities dj, a (possibly minimal) 
spanning set for functions in which to express multiple polylogarithms with singularities 
only in these Oj, provided such a set exists. The approach we present is rather generic and 
can be applied to any expression involving multiple polylogarithms. This is made possible 
by using results closely related to work of Goncharov, Spradlin, Vergu and Volovich [60], 
which in turn was inspired by the theory of (mixed Tate) motives 2 , and in particular by 
using a certain tensor calculus associated to iterated integrals, which is called "symbol 
calculus" in the following (the name "symbol" originating from [60] and from ref. [89]), 
and which we will review in the next section. 

An important remark is that the construction of a "symbol" seems to be a rather 
special case of a very general construction by Chen [22], where it appears as the image 
of an iterated integral as a 0-cocycle in the so-called "bar construction" attached to, say, 
X equal to the projective line minus a number of points (more generally, the construction 
has been investigated for a hyperplane configuration [90], §3), and it lands in the n-fold 
tensor product of the vector space of 1-forms on the underlying space X. Moreover, Chen 
characterised the image as the formal words in these 1-forms satisfying a natural integra- 
bility condition. Therefore, it would seem appropriate to reflect this in the notation for 
this object, e.g. as "Chen symbol". A polygon attached to an iterated integral enjoys the 
useful property that it gives a very concise way of explicitly producing integrable words, 
i.e. (Chen) symbols, of that kind. 

As an application, we restrict ourselves in section 6 to a specific subclass of multiple 
polylogarithms that are of particular importance in applications in high-energy physics. 
These so-called harmonic polylogarithms (HPL's) H(a; x) were first singled out and thor- 
oughly studied in ref. [25]. HPL's correspond to a special case of the iterated integral 
defined in eq. (2.1) where a\ € {—1,0, 1}. More precisely, they are defined via 

H(a;x) = (-l) k G(a;x), ^€{-1,0,1}, (2.14) 

where k is the number of elements in a equal to (+1). Many one and two-loop Feyn- 
man integrals can be expressed in terms of HPL's up to weight four and generalizations 

2 Let us point out that this is far from being the first exhibit of a direct connection between mixed Tate 
motives and mathematical physics, as such a relationship has been explored, e.g., by Kreimer in work with 
Bloch and Esnault [82, 83], such a connection was clearly apparent from letter correspondence between 
Broadhurst and Deligne [84] resulting e.g. in ref. [85], work of Belkale-Brosnan [86] or more recently by 
Brown [87] and others. One should also mention work of Connes and Marcolli [88] in this direction. 



- 5 - 



thereof [26, 27]. As harmonic polylogarithms are just a special case of the multiple poly- 
logarithms introduced at the beginning of this section, all HPL's through weight three can 
be expressed through classical polylogarithms. By contrast, similar to the general case of 
multiple polylogarithms, it is expected that HPL's of weight > 4 are no longer expressible 
in terms of classical ones alone. In section 6 we illustrate our technique by constructing a 
spanning set of harmonic polylogarithms in weight 4. 



3. Symbols and polygons 

The differential structure of multiple polylogarithms can be captured very well combi- 
natorially using a certain kind of decorated polygons with some additional structure, as 
developed in ref. [21], where they were called R-deco polygons. We note that there are 
related notions that had occurred previously in Goncharov's work, e.g. in refs. [8, 91]. 
There is an algebraic object attached to such a polygon, and hence to the corresponding 
multiple polylogarithm. This object, which has been dubbed a symbol in ref. [60], is an 
element in a tensor power of a certain vector space and contains a lot of information about 
the original function. 

3.1 An example in a nutshell 

In this subsection we give a quick idea of how, following ref. [21], one can associate to a 
multiple polylogarithm — or rather to an associated rooted decorated polygon — its symbol 
(we show in section 3.2 that this definition is equivalent to the definition given in ref. [60]). 
In the following subsection we then give a more detailed account of the construction. 

A multiple polylogarithm of weight n gives rise to a certain (n + l)-gon. As a "fore- 
shadowing" example, we first give the 4-gon P = P(c, b, a, x) attached to some weight 3 
multiple polylogarithm G(a, b, c; x) = —Li\^\(b/c, a/b, x/a) = 1(0; c, b, a; x): 



P(c,b,a,x)= i — > G(a,b,c;x) 



which comes equipped with decorations (in this order) c, b, a and x, the latter decoration x 
being for the distinguished root side (drawn by a double line in the picture) , and also carries 
information on the orientation of the polygon in the form of a fat vertex (which should be 
thought of as the "first" vertex, while the root side — adjacent to the first vertex — is the 
"last" side). 

In a first step, one lists all possible ways to draw the maximal number of non-intersecting 
arrows (an arrow is a directed line from a vertex of P to a non- adjacent side), which for 
an (n + l)-gon amounts to n — 1 such arrows, and one formally adds the resulting objects 
(the framing polygon being identical, but each equipped with a different maximal set of 
arrows) . 

In our example n = 3, such a maximal set contains n — 1 = 2 arrows, and there are 
precisely 12 different such sets, given by 
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X X X X 




b b b 



X X X X 




b b b 



X X X X 




b b b 



In a second step, to each such maximal set A of arrows in P, we associate a rooted 
tree (as the tree dual to the polygon dissection defined by the arrows) whose decorations 
are (decorated and rooted) 2-gons. As an example, to the 4-gons in the last column above 
we attach 





b 



Any linear order on the vertices of such a rooted tree which is compatible in the sense 
discussed below in section 3.2 with the partial order on it (only the middle tree above is 
non-linear hence allows more than one such linear order) now gives a term in the symbol 
S(P) attached to P. In practice, this means that every branching in a tree contributes to 
the symbol by the shuffle of (the vertices that appear on each of) its branches (see below 
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for a more detailed description). 

Third step: Each of the 2-gons B in one of the linear orders on the vertices now is 

mapped via a suitable map fi to a rational function in the original decorations of the 

polygons (in the example a natural target space would be the function field Q(a, b, c, x) of 

y 

rational functions in the variables a, b, c, x). More precisely, if B = where y 

X 

denotes the root decoration, then we map B to ^(B) = 1 — y/x, provided x / 0, and to 
Ijl(B) = y otherwise. 

Last step: Fixing the signs. We need to invoke a sign for the individual elementary 
tensors, and this sign is determined by using the number of backward arrows in a dissection. 
In order to see this quickly, it is convenient to "break up" the polygon at its first vertex 
(in the pictures it is typically indicated by a bullet). Then we "roll out" the sequence of 
sides and arrange it as a line from left to right, starting with the first vertex and ending 
with the root side; dissecting arrows inside the polygon will be stretched out (in a way that 
they still do not intersect). We give it for the third example above: 




Now a backward arrow is one which, in the rolled-out version of the polygon, has its 
end point to the left of its starting point (i.e. points from right to left, like (3 above), while 
a forward arrow has it to its right (i.e. points from left to right, like a above). 

Here is a more formal definition: There is a natural linear order on the sides e±, . . . , e n 
of an n-gon as above, starting with the non-root side e\ incident with the first vertex and 
ending with the root side e n (in the example above it is the linear order given by the 
sides ei, . . . , e4 decorated by c, b, a and x, and the vertices v\ = e n n e\ (the first vertex), 
V2 = e\ n e2, V3 = e2 H e3, V4 = e^ n e^. This induces a linear order on the vertices Vj 
which arise as the intersection of ej and ej-i (indices taken modulo n), where the first 
vertex is the smallest element in that order. Then a (non-trivial) arrow can be encoded 
by a pair (vj,ek) with k £ {j — 1, j}, and it is backward if and only if k < j — 1. With 
these notions, the sign attached to a polygon P with a maximal arrow set A is given by 
(— 1 )#{ backward arrows of A} ^ j n three examples discussed above in more detail we get two 
backward arrows for the first maximal dissection of the square and one backward arrow for 
the remaining two dissections. 

Putting all of the above ingredients together and writing ta for the tree dual to the 
maximal set of arrows A, and (ta, -<) for its partial order, the final formula for the symbol 



-8- 



S(P) of an n-gon P is 

g^pj ^ ^ ^ #{ backward arrows of A} 



E 



max sets A 
of arrows in P 



linear orders A 
compatible with (t^,-<) 





As an example, the first and third of the three maximal sets of arrows above give 




x 



and 



cab 

respectively, while the middle term (corresponding to a non-linear dual tree, i.e. a dual 
tree with branchings) contributes via the shuffle product of the two branches 

c ^ b a ' 

=-K*^) - ^K - ^^) 



__ (l _£) 8(l _«) 8(l _£)_ (l _ : 

where we introduced the symbol for the shuffle product 



aUIb = a®b + b® a . 



(3.2) 



Motivation and justification of this assignment has been given to an extent in ref. [21], 
where it forms part of an expression arising from the well-known bar construction in al- 
gebraic topology applied to a differential graded algebra on the polygons above (which in 
turn is motivated by certain algebraic cycles originally studied by Bloch [92] and Bloch- 
Kriz [19]). For an earlier appearance of a very similar structure (called the <g) m -invariant 
there), see ref. [8], §4.4. 

To summarize: an important part of the differential structure of a weight n multiple 
polylogarithm is captured by a certain decorated (n + l)-gon. More precisely, if the argu- 
ments of the multiple polylogarithm are expressed in terms of variables/constants xi, . . . , 
x m for some m, the polygon is an (n + l)-gon with decoration by simple expressions in 
xi, . . . ,x m ; now to this (rooted decorated oriented) polygon there is attached in a natural 
way an expression (its "symbol" ) in V® n where V is a finite rank submodule (it might be 
convenient for the reader to think of V as a finite dimensional vector space) of the space 
Q(xi, . . . , x m ) x (of infinite rank), i.e., the invertible rational functions in the variables 
X\i ■ ■ ■ i x m . 
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3.2 Rules of symbol calculus 

Roughly, a symbol is a formal sum of elementary n-fold tensors a± ■ ■ ■ a n , and one works 
in each tensor factor as with (a refined form of) dlog terms. In other words, each factor ai 
in a tensor product is tacitly understood as 

d l ogai = ^l. (3.3) 

Furthermore, we use shuffle products and the following rules (essentially boiling down to 
multilinearity, but in an unusual form, as we pass from multiplicative to additive notation): 

• Distributivity. 



C®(a-b)®D = C®a®D + C®b®D (3.4) 

and consequently 

C <g> a a <g> D = n(C <g> a <g> D) , neZ, (3.5) 

where C and D denote fixed elementary tensors. Note that n here is a coefficient 
rather than part of the first tensor factor; in particular, putting n = we see that 
C ®1® D = 0. 

• Neglecting torsion. 

We will "work up to torsion" , which means that we will put 

C®p n ®D = 0, neZ, (3.6) 

for p n an n-th root of unity. 

• Shuffle product. 

An important property of the symbol is that it preserves products: more precisely, 
it maps the product of two multiple polylogarithms to the shuffle product of their 
respective symbols, i.e. 

<S(G(ai, . . .,a r ;x)G(b 1 , . . .,b 3 ;y)) =S(G(a 1 , . . . ,a r ;x)) UIS(G(b 1 , . . .,b s ;y)) 

(3.7) 

where III is the symbol used for the shuffle product of two tensors, defined on 
elementary tensors by 

(«i ® . . . ® a ni ) III (a ni+ i ® . . . a ni+ „ 2 ) = ^ ® ••• ® a a -i( ni+n2 ) , 

<re£(ni,ri2) 

(3.8) 

where £(7)4,712) was defined in eq. (2.4). For more details on shuffle algebras, we 
refer to appendix A. We note that, on the left hand side of eq. (3.7), the shuffle 
permutations are applied to the arguments of the two functions (cf. e.g. eq. (2.3)), 
while on right hand side one shuffles the tensor factors instead, in a completely 
analogous fashion. 
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Note that eq. (3.7) is a rather non-trivial fact, as one can already see in the first 
non-obvious case: 

S(G(a; x)G{b; x)) = S(G(a, b; x) + G(b, a; x)) = 

= ((i->(i-f)-(i-><i-^)+(i-f)0O) 

+ ((i-f)«d-f) - (i-f)®d-|)) + (i-s)®(i-5)). 

which agrees with S(G{a; x))UIS(G(b; x)) = (l - § ) ® (1 - f ) + (l - f ) <g> (1 - f ), 
due to cancellations of terms. 

We will encounter expressions which involve both tensor and shuffle products — in 
order to avoid writing many parentheses, our convention is that a shuffle takes prece- 
dence over a tensor, i.e. 

a<g>6IHc = a<g> (6IIIc). (3.9) 
Furthermore, we abbreviate elementary tensors with the same factors as follows: 

q® n = a® • • • ®q . (3.10) 

n times 

• Refined dlog terms. 

We emphasise here already, though, that we will not treat d log c for a rational con- 
stant c as zero (as opposed to ref. [60]) since we would lose a lot of important in- 
formation this way. Instead we extend the above calculus to rational numbers in 
complete analogy with the above; so we have, e.g., 

C®2 m -3 n -x- 5 ®D = m(C®2®D) + n(C <g> 3 <g> D) - 5 (C x <g> D) . (3.11) 

• Root decoration annihilates: 

Since G{. . . ; 0) = 0, we also need to put S{G{. . . ; 0)) = 0, and this indicates that we 
can (and will) ignore polygons whose root side is decorated by 0. 

Linear orders of a tree. For a rooted tree T, which we view without a fixed 
embedding into the plane, hence e.g. we consider as equal the two trees 

and / \ 

V2 V2* \vi 

There is a natural partial order -< on its vertices Vj (j G J), given as follows: the root 
vertex vq ~< Vj for any j G J, and Vj -< Vk for vj / Vo if and only if there is a direct path 
from root to a leaf passing first through Vj and then through Vk ■ 

A linear order on the vertices of T which is compatible with the order -< is a sequence 
(vo,Vj 1 , ... ,Vj r ) of all the vertices Vj (j G J) such that Vj i ~< Vj k implies j% ^ 3k- (This 
means that if two vertices are in a relation with respect to the partial order, then they 
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should be related in any compatible linear order in the same way, while if they are not 
related in the partial order, there is no condition for how they should be related in that 
linear order.) In the example, there are precisely two linear orders which are compatible 
with the partial order, as the root vertex always comes first: (vo, v\, V2) and (^0,^2,^1)- 

Definition of a symbol. Now we are ready to give a complete definition of the symbol 
attached to a (rooted decorated oriented) (n + l)-gon P with decorations (ti, . . . ,t n ,x), 
and then extend it by linearity and shuffle product to any sum of (products of) polygons, 
hence also for multiple poly logarithms: 

S(P)= (-l)#i backwardarrowsof ^ ^ »{nt ,X ) ® ■■■®^n X ) , 

max sets A linear orders A 

of arrows in P compatible with (t~j±,-<) 

(3-12) 

where the 2-gons ir u ' are determined by the maximal dissection A together with the linear 
order A which is compatible with the partial order -< on ta, the dual tree of the dissection 
A, in the manner given above in the second step of section 3.1 (i.e. for each 2-gon arising 
from the dissection of A there is a vertex of ta decorated by that 2-gon, and for any two 
2-gons that are adjacent there is an edge in ta connecting the corresponding vertices. 

Integrability condition. A very useful property of the rooted decorated polygons, 
found by the second author in collaboration with F. Brown and A. Levin, is that each 
polygon (or rather its symbol) satisfies a certain integrability condition. Indeed, an arbi- 
trary sum of elementary tensors does not necessarily lie in the image of the symbol map. 
Instead, it was pointed out in ref. [90], making explicit in a special case the very general 
approach of Chen [22], that a necessary and sufficient condition for a symbol 

S= °i ® ■ ■ ■ ® co im (cj€Q), (3.13) 

I=(ii,...,i m ) 

to be integrable to a function is that 

^ cj uJi-t <g> ■ ■ ■ <g> (u) i:j Au^ +1 ) <g> ■ ■ ■ <g> uj im = for all 1 < ,7 < m, (3.14) 

I=(i 1 ,...,i m ) 

where Wj . A - +1 denotes the usual exterior product of two differential forms. We rewrite 
this for our purposes as 

^ cj [dlogw^. A dlogu ij+1 ] uj h (g) ■ ■ ■ (8)0;^. ® (g) ■ ■ • (g) uj im = 0, (3.15) 
J=(ii,...,i m ) 

where the hats indicate that we omit the corresponding factors in the tensor product. As 
an example, we indicate the statement for G(a,b;x), whose symbol 

5(G(a, 6 ;x)) = (l - f ) e (l - ?) - (l - f ) « (l - |) + (l - ?) « (l - ?) (3.10) 
satisfies 

dlog(l-|)Adlog(l-^) - dlog(l-|)Adlog(l-^)+dlog(l-^)Adlog(l-|) =0, 

(3.17) 
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where we recall that d\ogx = dx/x. Indeed, writing 

«\ dy — dj3 



d\og{\ 



the left-hand side of eq. (3.17) becomes 

dy - db x 



dy - db x ^dy - da * 
y =o y - a 



-b 



y=0 



y -b 



A 



dy — da 



(3.18) 



— a 



y=0 



+ 



dy — da 



y -a 



A 



dy — db 



y-b 



y=o 
(3.19) 



and we find, e.g., that the coefficient of dx A da is given by 



1 



-1 



1 



-1 



+ 



1 



1 



x — b x — a x — b b — a x — a a — b 
The coefficients of dx A db and db A da vanish in a similar way. 



= 0. 



(3.20) 



Relationship to the symbol of ref. [60]. In ref. [60], the Goncharov, Spradlin, 
Vergu and Volovich use the differential equation for multiple polylogarithms recursively to 
arrive at the definition of a symbol. More precisely, if F : C n — > C denotes a complex 
valued function depending on n complex variables Xk, 1 < k < n, the authors of ref. [60] 
define the "symbol of the transcendental function" F in the following recursive way: if the 
total differential of F can be expressed in the form 



dF = Y^ FidlogRi, 



(3.21) 



where Fi and Ri are functions of the variables x^, and Ri are moreover rational functions, 
then the symbol of F is defined recursively via 



S(F) = £ S(Fi) ® Ri ■ 



(3.22) 



In the case where F is a multiple polylogarithm, we can write down the differential of F in 
an explicit form. For example, in the special case where all the arguments of the multiple 
polylogarithm are generic (i.e., they are mutually different and do not take particular 
values), we obtain [2] 



dG(a n -!, . . . ,ai;a n ) = V"G(a n _i,.. . ,a;,. . . , a^, a n ) dlog I — — 

f-f \a-i-ai-i 
t=i 

The symbol of G(a\, . . . , a„_i; a n ) is then defined in the form 



(3.23) 



n-l 



S(G(a n -i, . . . , ai; a n )) = ^<S(G(a n _i, ai; a n )) 



i=i 



O-i — Qj+1 
O-i — O-i-l 



(3.24) 



The symbol we just obtained looks seemingly different from the definition we gave in 
eq. (3.12), which consists in summing over all possible maximal sets of arrows of the 
polygon P(a\, . . . , a n _i, a n ) associated to G(a n -i, . . . , a±; a n ). In the following we show 
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that the two definitions are equivalent up to a rearrangement of the terms in the sum, and 
hence give rise to the same tensor. 

Let us consider the n-gon P = P(a\, . . . ,a n ) (i.e. with sides decorated by ai, the 
last one a n decorating the root side). We will show that the symbol of P satisfies the 
recursion (3.24). For simplicity, we will concentrate here on the case of generic decorations. 
Let Ap be the set of all linear orders on the dual tree attached to any of the maximal sets 
of arrows of P. Then there is an obvious bijection between the terms in the double sum 
in eq. (3.12) and the elements in Ap. We can partition Ap by collecting all those linear 
orders into a subset which share the same last 2-gon that decorates the last vertex of this 
linear order. This partitions Ap into a priori 2n subsets, as those last vertices correspond 
precisely to the 2-gons that we can cut off from P. 

Note that cutting off the last 2-gon in a linear order on a maximal dissection corre- 
sponds to contracting the associated edge in the dual tree. Note also that, clearly, the last 
vertex must be a leaf of the (rooted) dual tree, and hence each last 2-gon necessarily cuts 
off two successive sides of P. 

Remark: For each n-gon P with n > 2, the three 2-gons involving the root side of P 
can never become the last one in any linear order in Ap . More explicitly, these are the 

two 2-gons K^^y and vj^^' . The former can arise only by cutting off the root side, 

a»-i ai 

while the latter can arise both by cutting off the root side and by cutting off the first side 
and the first vertex. As a consequence, Ap partitions into only 2n — 3 non-empty subsets 
(of same cardinality) of the above type. 

In view of the above, it is clear that any such subset indexes exactly the linear orders 
Ap on the (dual trees of the maximal dissections of the) subpolygon P of P which is 
obtained from P by cutting off a fixed 2-gon, followed by contracting the dissecting arrow 
to a point. 

There are typically two ways of cutting off a 2-gon in which such a subpolygon P can 
occur: cutting off the 2-gons (« = 2, . . . , n— 1) leaves complementary subpolygons 

Pf~ which are identified with Pj upon contraction of the dissecting arrow. Note that the 
two 2-gons will give terms of opposite parity, as precisely one of them will be a forward 
arrow. The only exception arises from cutting off *^=y , which corresponds to a forward 

ai 

arrow, for which only one complementary subpolygon can occur. 

In summary, we get: 
Claim: There is a bijection of sets 

Ap ^ y L pt x {^]\ u y L p _ x {^]\ . (3.25) 
i=i y ai J i=2 y ^ J 

Moreover, the sign of a maximal set of arrows of hp agrees with the sign of the correspond- 
ing maximal set of arrows in A p + and is opposite to the sign of the corresponding maximal 
set of arrows of Ap-. 
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All we need to note here is that the above remark ensures that both A„± and A5- 
are empty, so are left out at the right hand side, and that cutting off a 2-gon of the form 
corres P on ds to a forward arrow, hence contributes a sign +1 to the maximal 

a i — l 

dissection, while corresponds to a backward arrow, hence contributes a sign — 1. 

a. 

Due to the bijection (3.25), we can rewrite eq. (3.12) by first summing over all sub- 
polygons P^, followed by a sum over all possible elements in As±. The inner sum then 

i 

evaluates to the symbol of the subpolygon Pf 1 , and we are left with 

rt— 1 n— 1 

i=l ai i=2 ai 

where the relative minus sign arises because, as discussed above, Ap± contribute with 
opposite signs. After identification of and P~ , eq. (3.26) agrees with the recursion (3.24) 
modulo the additivity of the symbol. In order to finish the proof, we need to show that also 
the bases of the recursions are the same. It is indeed easy to check by explicit computation 
that, e.g., the symbol of G(a,2, 0,1', 0,3) obtained from the recursive definition (3.24) agrees 
with the symbol obtained from our polygon construction, eq. (3.12). 

We note that in ref. [60] the d\og of a constant is put to 0. Although this seems rather 
natural (and turns out to be sufficient in several cases), we advocate to use a refined version 
of this (which is what is typically used when working with symbols in a number field) : for 
each element of a set of multiplicatively independent elements in a given number field one 
can choose a logarithm independently but then the logarithm of any product formed from 
those elements is determined. For example, we will see in section 6 that in the context 
of harmonic polylogarithms the only constants that need to be treated in this fashion are 
powers of 2, and hence it is sufficient to think of "2" as an irreducible element. The reason 
for considering this refined version is that it is very helpful for recognising functions from 
which a given symbol might originate. In particular, it has proved to be very useful, e.g., 
when "recognizing" HPL's for keeping track of terms which come from a (shuffle) product 
of polylogarithms, see section 6. 

While the symbol of a multiple polylogarithms obtained by considering the maximal 
set of arrows of the associated decorated polygon is equivalent to the symbol obtained 
from the recursive definition (3.24), we believe that both approaches have their virtues. 
While the latter might be easier to implement into a computer program in general, it is 
strictly speaking only valid in the case of generic arguments of the polylogarithms. Indeed, 
if the arguments are non-generic, we obtain divergences in the right-hand side of eq. (3.23), 
e.g. when cij = ctj+i for some i. It is then in principle necessary to resort to a careful 
regularization to deal with the degenerate cases [2] . The definition of the symbol based on 
the decorated polygons, being combinatorial in nature, avoids this problem by construction 
and allows to identify very easily the degenerate cases (they correspond, e.g., to arrows 
ending on a side whose decoration is zero), and to discard them from the start, avoiding 
in this way the need of regularization. Furthermore, as we will see in the next section, 
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the polygon approach has naturally built in the refined l d log- prescription', because the 
combinatorial nature of the construction does not make a distinction between constants as 
e.g. "2" for which one might be typically tempted to define cHog2 = 0. 

Symbols for classical polylogarithms. The polygons attached to classical polylog- 
arithms Li rn (x) = —G( 0, . . . , , 1; x), are given by decorations x (for the first side) and 

m— 1 terms 

(for the remaining non-root sides) as well as 1 (for the root side). Their attached symbol 
consists of (the negative of) a single elementary tensor, in fact we have 

S(Li m (x)) =- ((l-x)® x®-y®x ) , (3.27) 

m— 1 factors 

where we have m factors ("weight" m) on the right hand side. (Note the parentheses which 
separate the coefficient, here —1, from the actual tensor, to avoid a misinterpretation as 
(x — 1) (8 x (8 ■ ■ ■ <8 x.) 

Such tensors have long been considered in connection with functional equations of 
polylogarithms — in fact, Zagier [5, 58] has given a criterion for such equations built on 
those tensors, which has been used (cf. ref. [59]) to find the first non-trivial equations for 
Lis and Lij (beyond weight 7 none are known), and the corresponding expressions for 
multiple polylogarithms are important already in Goncharov's early work (e.g. [7]) where 
he generalises the underlying tensor algebra considerably. 

4. A simple example 

The symbol attached to G(—l,l;x). In this section we illustrate the fact that the 
symbol calculus provides a convenient tool to detect functional equations among (multiple) 
polylogarithms, on the example of G(— 1, l;x) (which happens to coincide with the HPL 
—H(—l,l;x)). Even though we could of course immediately apply eq. (2.13) to express 
G(— 1, 1; x) in terms of classical polylogarithms, we will derive a similar functional equation 
using the tensor calculus introduced in the previous section. The multiple polylogarithm 
G(— 1, 1; x) is associated to a trigon, 



X 



G(-l,l;x) <— ► P(l,-l,x) = 




The dissection of the trigon can easily be translated into the tensor associated to the 
polylogarithm, 




®K^) -K^M 
-ii 1-11-1 
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The last line allows to read off the symbol of G(— 1, 1; x), 



S(G(-l,l;x)) 




= (1 + x) <g> 2 + (1 - x) <g> (1 + x) - (1 - x) <g> 2 . 



Before turning to the question of how to attach a function to this symbol, let us 
briefly comment on how this symbol could have been obtained by using the recursive 
definition (3.24). Using eq. (3.23), we obtain 

dG(-\,l\x) =G(-l;z)dlog(l-(-l)) + G(l;x)dlog 

= G(-l;x)dlog2 + G(l;x)dlog(l + x) - G(l;x)dlog2. 

The three terms in the last line of this equation are in one-to-one correspondence with 
the three terms in the symbol in eq. (4.2). Note, however, that we had to treat all the 
arguments of the (three- variable) function G(», •; •) as generic, and to use the refined l d log- 
prescription', i.e. a! log 2 ^ 0. Putting to zero all the d log 2 terms is equivalent to putting 
to zero all elementary tensors in the symbol where a factor inside a tensor product is 
constant [60]. As we will see below, we prefer to keep these terms as they provide us with 
valuable information about the function that should be associated to the symbol. 

Attaching a function to a symbol. Since every multiple polylogarithm of weight two 
can be expressed as a combination of classical polylogarithms, we make the ansatz that 
G(— 1, 1; x) can be written, up to an additive constant, in the form 

^CjLi 2 (/;(x)) + J2 c i k l °g(9j(x)) log(Mx)), (4.4) 

i j,k 

such that the tensor associated to this expression corresponds to the tensor in eq. (4.2), 
where c, and Cjk are rational numbers and /«, gj and hk £ Q(x) are rational functions. We 
subdivide this problem into smaller ones by postulating that we can distinguish between the 
three different contributions in eq. (4.2). By using a procedure suggested in ref. [60, 75] 
we can distinguish the first sum from the second by projecting the respective symbols 
onto their symmetric or alternating parts: each term in the second sum will give zero 
contribution for the latter one, while each summand in the first sum will give a non-zero 
contribution. Indeed the tensor associated to a product of logarithms is totally symmetric, 
and hence its contribution to the tensor vanishes if we project onto the antisymmetric 
component of the tensor in eq. (4.2). 

Preparatory steps: decomposing tensors into symmetric and antisymmetric 
parts. We recall that, for a vector space V (over C, say, or more generally over a field of 
characteristic ^ 2) there is a direct sum decomposition V <8> V = {V QV) ® {V f\V) (other 
notations, as used e.g. in refs. [5] or [7], are V V = Sym 2 (V), and V A V = f\ 2 V), and 
V V is generated by a b (for some a, b G V), while V A V is generated by a A b where 



V(-i)-V 
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we introduce the following rather standard notations for symmetric and antisymmetric 
tensors, 

aQ b = a<g>b J t-b<g>a, 

(4.5) 

a Ab = a 6 — 6 a . 

Back to the example. With this notation, the decomposition of a generic elementary 
rank two tensor (i.e., a b for some a, b £ V) into its symmetric and antisymmetric 
components can be expressed as 

a<g>b= ^ (a 6) + ^ (a A 6). (4.6) 

Concentrating solely on the antisymmetric component of eq. (4.2), and using the antisym- 
metry of the wedge product, which, e.g., induces 2 A 2 = 0, we obtain 

(l + x)A2 - (l-x)A2 + (1 - x) A (1 + x) 

1-x l + x ( l + x\ l + x (4.7) 
= — — A — — =1 — A — — . 



As the tensor associated to a product of logarithms does not have an antisymmetric com- 
ponent, eq. (4.7) suggests that it is the antisymmetric part of the tensor associated to some 
sum of dilogarithms, and from eq. (3.27), it is easily identified as the antisymmetric part 
of<S(-Li 2 (^)). 

Having identified the dilogarithm contributions to G(—l,l;x), we can proceed via a 
bootstrap procedure and subtract off this contribution, leaving only a totally symmetric 
tensor 

S (g(-1, 1; x) + Li 2 (±±*) ) = 2 (1 + x) - \ (2 2) 

= S flog 2 log(l + x) - - log 2 2 

Fixing the constant. We have shown that the tensor associated to G(—l, l;x) equals 
the tensor associated to the combination — Li 2 (^-p) +log2 log(l + x) — \ log 2 2. It would 
be premature, however, to conclude that both expressions are equal, but they are equal 
only up to an additive constant independent of x. Indeed, specializing to x = 0, and using 
the fact that G(-l, 1; x = 0) = and Li 2 (l/2) = - \ log 2 2, we see that 



|G(-l,l;x) 
Thus, we obtain 



-Li 2 ( j + log 2 log(l + x)-- log 2 2 



_ 7T 2 

i " 12 

x=0 



(4.9) 



G(-l, 1; x) = -Li 2 ( + i g 2 log(l + x) - \ log 2 2 + ^ . (4.10) 



2 J ° ' 2 ° 12 

Note that this additive constant is not detected by the symbol, because 

S(t?) = _5(log 2 (-l)) = . (4.11 
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5. Integrating symbols: an algorithmic approach 

In the previous section we illustrated how the symbol calculus can be used to derive a 
functional equation among poly logarithms. While in that weight two example all the steps 
were easily carried out by hand, an algorithmic approach is desirable in more complicated 
cases. In this section, we present our approach that often allows to integrate a symbol of 
a (transcendental) function, i.e., to find a function F, written as a linear combination of 
(products of) multiple poly logarithms, whose symbol matches a given tensor S, which in 
the following we always assume to satisfy the integrability condition (3.15). 

From the example of the previous section it should be clear that the main challenges 
to address when going to higher weight w are 

1. choosing appropriate arguments of the desired function types (as a few examples of 
function types of weight four, we list L14, Li2,2 or Li2 log log) such that their symbols 
span the vector space in which the tensor S lies; 

2. finding the generalization of the decomposition of weight two tensors into symmetric 
and antisymmetric parts (indicated in the simple example of weight two in section 4) 
to higher weights. This problem was addressed up to weight four in refs. [60, 75]. 

Let us assume that we have a linear combination S (with rational coefficients) of 
elementary tensors where the factors in each elementary tensor are rational functions of 
some variables x±, . . . , x n . In the following we assume the tensor to be of "pure weight" w, 
i.e., each elementary tensor is assumed to have the same number of factors. Without loss 
of generality we can then assume that S takes the form (all sums are assumed to be finite) 



where Ri e (xi, . . . , x n ) (1 < it < m(S) for some m(S) determined by the initial tensor S) are 
rational functions in the variables X{ and the Cjj ^ are rational numbers. Distributivity 
(cf. eqs. (3.4, 3.5)) then implies that, without loss of generality, S can be rewritten (with 
new constants Cj lt ___j w 6 Q) as 



where ttj = ttj(xi, . . . , x n ) (1 < j < K for some K) are suitably chosen rational functions 
which are multiplicatively independent (i.e., there is no non-trivial relation of the form 
IT^=i T^j = il> f° r r j In practice, we achieve this by simply factorizing the rational 

functions Ri(xi, . . . , x n ) in eq. (5.1) into irreducible polynomials over Q, say (i.e., polynomi- 
als in Q[xi, . . . , x n ] that cannot be written as the product of two non-constant polynomials 
in Q[xi, . . . , x n ]). Let us denote the set of these polynomials by Ps = {tti, ■ ■ ■ , kk}, which 
will constitute our building blocks in the following. The symbol can then be seen as an 
element of the tensor algebra over the Q- vector space generated by the formal basis vectors 
in the set Ps (more precisely, we should consider it as an element of the wth grading of 




(5.1) 




(5.2) 



Weight 


Basic function types of pure weight 


1 


log X 


2 


Li 2 (x) 


3 


Li 3 (x) 


4 


Li 4 (x), Li 2j2 (x,y) 


5 


Li 5 (x), Li 2;3 (x,y) 


6 


Li 6 (x), Li 2j4 (x,y), Li 3)3 (x,y), Li 2)2j2 (x, y, z) 



Table 1: Indecomposable multiple polylogarithms of pure weight < 6. 



the tensor algebra over the Z-module (± IX;=i I r j ^ ^))- ^ur goal is now to find a 
function, say F, that is a rational linear combination of (multiple) polylogarithms (and 
products thereof) whose arguments are rational functions in the x« such that S(F) = S. 
The procedure to achieve this proceeds in two steps: first we have to decide on the types of 
functions that should appear in F, and then we have to concoct suitable rational functions 
in the guments of these functions such that for some linear combination of these 

functions the resulting expression fulfills the condition S(F) = S. Note that this latter 
step is not algorithmic in general, as it may involve finding arguments for the functions 
that have singularities outside Ps- 

5.1 Choosing the types of functions 

Our first goal is to construct a set of function types (our 'basic types') out of which we 
can construct our candidate function F. This involves two steps, and we want both the 
functions and their arguments to be 'as simple as possible', but we need to take into 
account that all the possible function types one can write down for a given weight are 
related by an abundance of functional equations. The main criterion we will use in the 
following is that a function type that can be written as a product of lower weight function 
types is 'simpler' than a function type of pure weight/transcendentality (i.e., a function 
that cannot be written as a sum of terms, each of which being a product of lower weights). 
Furthermore, we are guided by the conjecture (which the second author learned many 
years ago from Goncharov) that a multiple polylogarithm Li mi) ... jmjfc with rrij = 1 for some 
j can be expressed in terms of multiple polylogarithms where no index is equal to 1. This 
conjecture suggests to put a restriction on the function types of pure transcendentality 
that can appear for a given weight. Furthermore, the shuffle and stuffle relations provide 
us with further constraints. As an example, we can deduce from 

Limi,m2 V) Lira 2 ,mi (y, x) = Li mi (x) Li m2 (y) - Li mi+m2 (x y) . (5.3) 

that we can hence ignore Li m2jmi with ra 2 < mi. For low weights, the corresponding sets 
of (presumably independent) functions which are indecomposable, i.e. cannot be written 
in terms of products of lower order functions, are given in Table 1. 
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5.2 Finding the arguments 

Having at hand a suitable set to construct the basic function types from, we still need 
to find the arguments of these function types. In the context of particle physics it has 
proved helpful to use guidance from educated guesses, motivated by physical constraints 
(cf. refs. [64, 65]), to construct the symbol and/or the functions expressing the desired 
physical quantities. To see how one might proceed even without any such guidance, let us 
concentrate first on classical polylogarithms only. 
We start by defining, for P$ = as above, 

P s = P s U P! s , (5.4) 

where P' s is the set of all prime factors that appear in 7Tj ± TTj and 1 ± 7Tj, Wj, Ttj £ Ps- 
Let us denote the elements of Ps by 7fj. Since S is constructed out of the irreducible 
polynomials 7Tj £ Pg C Ps, it is perhaps natural to hope that all arguments appearing in 
the polylogarithmic expressions needed for S can be written in the form 

R t 1 ,...,n k ( x U---i x n) = iTT^Oci, • • • ,x n ) ... 7f£ fe (xi,...,x n ), (5.5) 

where ni, . . . , are integers. Let us denote the set of these functions nk ( x i> ■ ■ • > x n) 

by IZs, i.e. this is, up to sign, the multiplicative span of the 7tj 6 Ps- Note that in practice 
it is often enough to consider IZs to be the span of only a subset of the polynomials in 
P. Finally, note that the set TZs carries a group structure, given by the multiplication of 
rational functions. 

Choosing arguments for classical polylogarithms. However, not all of these func- 
tions are good candidates for arguments of polylogarithms. Indeed, if for example such 
a function appears as an argument of a classical polylogarithm, then by eq. (3.27) for 
R = Rn lt ...,n k ( x ii • • • > x «) we can write- 

S(Li n (R)) = -(l-R)®R®---®R. 

v -v y (5.6J 

(n— 1) times 

It is now easy to see that if we want this tensor to be an element of the tensor algebra of 
the vector space generated by the set Ps, then we need to impose the constraint 

1-ReKs- (5.7) 

Let us introduce the notation 

= {R G TZ S | 1 - R G TZs] C 1Z S - (5.8) 

It follows that, for R € T^g\ the symbol of Li n (i?) (n > 1) is a linear combination of 
tensors of the form 7f^ 1 <g> . . . <g> 7f^ n . Hence all the rational functions in the set TZg 1 ^ are good 
candidates for arguments of the classical polylogarithms that can appear in our function F. 
Note that 1Z^ is no longer a group in general. It can however be given some more structure 
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by considering the following action of the symmetric group S3 on rational functions, defined 
for a rational function R and rational functions Oi of one variable by 



a 1 (R) = R, a 2 (R) = l-R, a 3 (R) = l/R, 
(j 4 (R) = 1/(1 - R) , a 5 (R) = 1- 1/R, a e (R) = R/{R- 1) . 



(5.9) 



Note that the aj form a group (under composition of functions) isomorphic to the permu- 
tation group S3 on three letters. It is easy to check that TZ^ is closed under this action of 
S3. As S3 is generated by the two elements 02 and (73, it is enough to check that is 
closed under these two maps. Closure under 02 is trivial by definition of T^g ■ To see that 
it is also closed under 0-3 we have to check that Vi? G 1 — C3(i?) G Us- Indeed, we 

have 

1 - a 3 (R) = 1 - l/R = -(1 - R) R' 1 G lis , (5.10) 
because of the group structure of Us- 

Choosing arguments for poly logarithms of depth > 1. So far we have only consid- 
ered classical polylogarithms, but in general we should also be able to make a sensible ansatz 
for the arguments of multiple polylogarithms of depths greater than one. In the following, 
we find it more convenient to work with the functions G 



mi,...,m k 



than with the functions Li 



mi,...,m k - 



defined in eq. (2.10) rather 
As the two function types are related via eq. (2.10), 
one can easily convert from one representation to the other. 

Let us consider a multiple polylogarithm of depth two, say C?2,2- We are hence looking 
for a pair of rational functions (i?i,i?2) G Us x T^S that are good candidates for the 
arguments of (?2,2- The symbol of G2,2(-Ri, R2) is given by 



S(G±2(R ■ rh:) - - [ I" ^ ) <V HI 



1 



1 - 

2 ^1 

f 



±-)®R 2 m 



1 

~ R2 
1 

R~2 



R2 

Ri 



R2 
Ri 



R2J R2 
l 

~R~2 



<g> II- 
>R 2 m 



R2 

' Ri 
1 

Ri 



® Ri III R 2 

(g) Ri <g> Ri 



1 



1 

R~i 



Ri 



recalling our notation for the shuffle products (see eq. (3.9)), 

A® B <g) CUID = A®B(g)C(g)D + A(g)B®D®C, 
A ® BUI (C ® D) = A ® B ® C ® D + A ® C ® B ® D + A ® C ® D ® B . 



(5.11) 



(5.12) 



Using the same reasoning as for classical polylogarithms, we see that the candidate argu- 
ments for multiple polylogarithms of depth two are pairs of rational functions from the 
set 



n 



(2) 



{(Ri,R 2 ) g n { s ] x n K s J \Ri - R2 g n s } 



(5.13) 
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An important consequence is that no new rational functions are needed to construct the set 
T^g \ but all the information is already contained in IZ^ The new set TZ^ then consists 
of pairs of elements of TZ^\ subject to the additional constraint that their difference must 
again be factorizable in terms of the same prime elements. Moreover, we saw that TZ^p is 
endowed with a natural action of the group S3, defined in eq. (5.9). It is hence natural to 

(2) 

ask for non-trivial symmetry groups that leave the set 1Z S invariant. First, it is easy to 

(2) 

see that the defining condition for TZ S is invariant (up to an overall sign) under swapping 
the two entries of any given pair. Second, the action of S3 defined in eq. (5.9) induces a 
(simultaneous on both factors) action on 1Z^ x 1Z^\ defined for <7j G S3 by 

{R 1 ,R 2 )^h{a i {R 1 ),a i (R 2 )). (5.14) 
(2) 

It is now easy to check that 1Z S is closed under this action. To see this, it is enough to 

(2) 

check that (Ji{Ri) — (Ti{R 2 ) £ TZs for i = 1,2 and whenever (Ri,R 2 ) G IZg '. Indeed, we 
have 



<72(Ri) - cr 2 (R2) = (l - Ri) - (l - R2) = -(Ri - R2) € TZ S , 
<ra(Ri) ~ ^(#2) = l/Ri - I/R2 = -(Ri - R 2 )Ri 1 R 2 1 e 1Z S , 



(5.15) 



where we used the fact that R\—R 2 6 1Z 5 and that TZs is a multiplicative group. Combining 
this S3 symmetry with the invariance under an exchange of arguments, here R\ and R 2 , 
we see that TZ^ is closed under the action of the group S3 x S 2 , defined for (a, p) G S3 x S 2 
by 

(R 1 ,R 2 ) H } (a(R p(1) ),a(R p{2) )) , (5.16) 

i.e. the factor S 2 simply acts as a permutation of the entries. 

The previous discussions for depths one and two readily generalize to higher depth. 
Our candidate arguments for the multiple polylogarithms of depth k are fc-tuples of rational 
functions from the set 

TZ ( g ] = {(#1, . . . , R k ) G t4 1} x . . . x TZ ( g ] \Ri-Rj£7Z s ,l<i<j <k}, (5.17) 

(k) 

and using exactly the same argument as in the depth two case, we see that 1Z S can be 
equipped with an action of the group S3 x S^, acting on (R±, . . . , R^) 6 TZ^ via 

(R 1 ,...,R k ) ( ^4{a(R p{l) ),...,a(R p{k) )) , (5.18) 
i.e. the factor S k simply acts as a permutation of the entries. 
5.3 Integrating the symbol (1) 

We now turn to the problem of integrating the tensor S that satisfies the integrability 
condition (3.15). In practice, such tensors could come from computing a Feynman integral 
in terms of multiple polylogarithms, or by computing its symbol by other means [63, 64, 
65, 67, 70, 73]. Our goal is to find a function F, more precisely a linear combination of 
(multiple) polylogarithms, such that S(F) = S. The considerations of the previous section 
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suggest writing down an ansatz for the set <3> of functions (i.e., function types, together 
with rational functions as arguments) to express F in. We assume the "folklore conjecture" 
that any functional equation of MPL's involving different weights should actually split into 
functional equations of pure weight, so we can assume S to be of pure weight. It is 
convenient to partition into subsets according to the weight w, say = {b\ w ^}i. 
Then we want to find rational coefficients c^...^ such that 

s + £ Cni2 s ( b ™ &£»>) + ...+ £ Cil .., w s (&£> . . . £>) 



11,12 
1U1+1D2 =11) 



(5.19) 



where we can assume that Wj > itfj+i for all z. In view of property (3.7), each term on the 
right hand side can be written as a shuffle of terms S (b^) . All the terms on both sides of 



this equation are elements of the grade w part of the tensor algebra over the vector space 
spanned by basis vectors labeled by the elements in Ps- We know from linear algebra 
that a straightforward basis for the tensor space of weight w tensors over the vector space 
spanned by Ps is given by {n^ (8). . . <8>7fj ro |7f^ £ Ps}- At this stage we have hence mapped 
the problem of finding a function F satisfying S(F) = S into a problem of linear algebra, 
more precisely the problem of finding the coefficients c% x ...i l such that eq. (5.19) is true. As 
we know a basis of the tensor algebra, we can just compute the coefficients by extracting 
and comparing the coefficient of each basis vector on either side of the equation and solve 
the ensuing linear system. Note that the function F obtained in this way is not unique. 
Indeed, the map S is non-injective and hence it would be possible to find a different function 
F' such that F — F' in the kernel of S. This issue will be addressed in section 5.6. 

5.4 A set of projectors 

Even though we have solved the problem of integrating the symbol in principle and we 
have reduced it to a linear algebra problem, the linear system one has to solve can be quite 
large. It is therefore preferable to break it down into smaller problems, for example by 
introducing a suitable filtration on the target space 3 which allows to successively solve the 
problem for the filtration pieces. Such a filtration would allow to separate the different 
contributions in eq. (5.19): to get the ball rolling, we would like to solve for the coefficients 
of the functions b\ w ^ without having to care about the product terms. This can be achieved 
by introducing a projector that sends to zero exactly all the product terms. 

Definition 1. Let V be a vector space. We define linear operators U w acting on elementary 
tensors of lengths w > 1 by Hi = id and for w > 2 by 

K w (ai <8> . . . <8> a w ) 

w - 1 (5-20) 
= — [II w _i(ai <8> ... <8> a w -i) <g> a w — n w _i(a 2 <8> • • • <8> a w ) <g> a±] . 



3 A filtration of a vector space V is a sequence {Vi}i<i<„ of subspaces of V such that Vi C V2 C . . . C 

V n = V. 
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Similar sets of operators acting on shuffle algebras have already been defined in refs. [78, 
93], differing from our definition only by an overall normalization (F. Brown, in a text 
in preparation on the representation theory of poly logarithms, uses essentially the same 
operators, referring to ref. [94], II §2, Prop. 7a). More precisely, the operators in refs. [78, 93] 
are defined by pi = id, and for w > 2 by 

(ai ® . . . ® a w) — Pw—1 (ai ® . . . <g> a w— l) ® 0"w Pw—1 (a 2 <8> . • • <8> a w ) <8> ai . (5.21) 

The exact correspondence between H w and p w is simply, for all w > 1, 

n w = - Pw . (5.22) 

It follows from refs. [78, 93] that the kernel of p w corresponds exactly to the ideal 4 generated 
by all shuffle products. Since II m and p w only differ by an overall normalization, we 
immediately arrive at the following 

Proposition 1. The kernel ofH w equals the ideal generated by all shuffle products, i.e., 
for every element £ in a shuffle algebra, II W (£) = if and only if £ can be written as a 
linear combination of shuffle products. 

In other words, the projectors Tl w are by construction such that they annihilate pre- 
cisely all shuffles. Conversely, if H w applied to some tensor £ does not vanish, then it is 
not possible to express £ as a linear combination of shuffle products. 

The reason for the normalization factor in Definition 1 is that it makes H w idempotent, 
in other words, U w is a projector. 

Proposition 2. For any w £ N , U w is idempotent, and hence a projector, i.e. 

U 2 W = U W . (5.23) 

Proof. In ref. [78], Lemma 1.2, and also in eq. (2) in ref. [93], it was shown that p w satisfies 
the identity 

w— 1 

w (ai <g> . . . <g> a w ) = (ai ® . . . ® a k ) III p w -k{ak+\ ® • • ■ ® a w ) . (5.24) 

fc=0 

We now act on both sides of this equation with p w . Since p w annihilates all shuffles, all 
the terms on the right-hand side vanish, except for k = 0. Hence we obtain 

w p w {ai <g> . . . <g> a w ) = p^ifli <g> . . . <g> a w ) , (5.25) 

and so 11^ = H w , after dividing both sides by w 2 . □ 

4 An ideal in a commutative algebra A is an additive subgroup I such that Va G A and Vb £ I, one has 
ab e 1. 
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5.5 Integrating the symbol (2) 

The projectors defined in the previous section can be used to improve the chances of suc- 
cessful integration of the symbol described in section 5.3. As the kernel of Tl w corresponds 
to all possible linear combinations of shuffles, its effect on eq. (5.19) is to remove all product 
terms, i.e., 

n w s = Y,an w s(b ( r ) ) , (5.26) 

i 

and so we can solve for fewer coefficients Cj without having to worry about the product 
terms. 

Having found the coefficients Cj, we are left with the determination of the coefficients 
of the product terms. In order to proceed via induction, we determine the behaviour of 
shuffles under tensor products of projectors (which are themselves projectors). We first 
study the instructive case of a projector Ug 8> Ibw which we apply to a shuffle 

(ai®---<g>afc)IH(&i<g>---<g>6fc/), (5.27) 

where k + k' = £ + £' = w (still w denoting the weight), with k > k! > 1, i > £' > 1 and 
k 7^ i. We can rewrite the shuffle in a form more suitable to applying IIg(g>n^/ by regrouping 
the sum (5.27) using deconcatenation (i.e. splitting into a left hand part and a right hand 
part) of the sequence (1, . . . , k) into two blocks (1, . . . , k\) and (k\ + 1, . . . , k\ + k 2 ), giving 

min(fc/) 

^ ((ai<g>---(g>a fe JIII(&i<g>---<g>^_ fcl )) (g> ((a kl +i<8)---<8)a k )UI(be-k 1 +i<8)---<8)b k >)) . 
fci=i 

(5.28) 

Applying lig <8 Tle> to this sum will annihilate each of the mm(k, £) summands individually, 
in fact it will annihilate the left hand factor except (possibly) for k± = £, since in the other 
cases there is a proper shuffle on the left which will be mapped to zero by He already, by 
Proposition 1. If the remaining case occurs, the corresponding right hand factor is mapped 
to zero by instead. For example, if (ki,k 2 ) = (3,3) and (£\,£ 2 ) = (4,2), we rewrite 

(ai ®o 2 a 3 ) III (bi <8 62 <8 63) = («i bH (h b 2 63)) ® («2 <8) a 3 ) 

+ ((ai®a2)in(6i®62)) ® (asIH^s) (5-29) 
+ ((ai a 2 ® a 3 ) HI (61)) (6 2 ® ^3) , 

where already the factor II4 of II4 <8> II2 annihilates the left part of each of the terms. 

More generally, for a partition A = (Ai, . . . , A r ) of w, we call the X-shuffte of a w-fold 
tensor a\ <8> • ■ ■ <8> a w the product 

(ai (8) • • • (8) a Al ) HI (a Al +i <8> ■ ■ ■ <8> aA!+A 2 ) HI ... HI (a Al+ ...+A r _ 1 +i (8 ■ • • <8> a w ) (5.30) 

and given any partition A' = (A' l5 . . . , A' s ) with s < r, the projector H\i will vanish on a 
A-shuffle, as we can again regroup its terms according to A' and find for each so combined 
summand at least one projector factor Uy which vanishes on the corresponding part. This 
sketches a proof of 
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Proposition 3. For two different non-increasing partitions X and X' of w, of length £(X) 
and£(\'), respectively, we have that a X-shuffte vanishes under Hy whenever £(X) > £{X'). 

This proposition suggests to define a sequence of subspaces, each of which contains 
the next one (hence forming a descending filtration). For this, we consider the standard 
lexicographic order y on non-increasing partitions of w given by (denoting [x\ the largest 
integer < x) 

H y (w - 1, 1) y (w - 2, 2) y . . . y (w - l|j , l| J) >- (w - 2, i, i) y 

y---y (2, (^^1 ) . 

w— 2 slots w slots 

As an example, in weight five we have the following non-decreasing partitions, ordered as 

(5) y (4, 1) y (3, 2) y (3, 1, 1) y (2, 2, 1) y (2, 1, 1, 1) y (1, 1, 1, 1, 1) . (5.32) 

From these, we form the descending filtration alluded to above by taking the span of 
all shuffles of partition types greater than a given type A, say, as 

T x = (A' - shuffles in V for all A y A', A + X 1 ) . (5.33) 

Hence is the zero space, J r (2,i,...,i) is generated by all (1, . . . , l)-shuffles (i.e. fully 

symmetric tensors), ^"(3,1,. is generated by all (1, . . . , l)-shuffles and (2, 1, ... , reshuffles 
etc. The proposition now guarantees that F\ C kerllA. Similarly, we put, for convenience, 
the "shifted" sequence 

JT^ = (A' - shuffles in V for all A t_ A') , (5.34) 

so that Jy + = T\ where A + denotes the immediate successor of A in the lexicographic 
order above. 

We proceed by induction on A with respect to the order starting with the shuffle 
A = (w). By the analysis above, we can put Si w \ = ^ i CiS(b\ w ^ , provided we have been 
able to solve for the q. This gives us the basis step for the induction. Now assume we 
have found an (integrable) tensor S\ "approximating" S in the sense that S — S\ G F\. 
Then we try to construct a "better approximation" S x + for the successor partition A + 
of A by finding integrable tensors T A + G Tr + = T\ such that S x + = S\ + T A + satisfies 
S — S x + G J~\+. If we are successful in finding such a , this finishes the induction step. 

We expect to be able to find such a 7\+ using (certain sums of products of) multiple 
poly logarithms, by selecting for each part A+ of A + an "indecomposable function type" 
of weight A+ (see table 1) and taking their product. In other words, we assume that the 
tensor T A + can be written, as a linear combination of A + -shuffles, in the form 

T A+ = £ c il ,..., il s(b^ ) )ui...Ws(bf } ), (5.35) 

h, —,k 



(5.31) 



with I = £(X + ) and A + = (Xf , . . . , A+). In the weight 5 example of table 1, the respective 
polylogarithmic functions we consider are 

(5) oLi 5 (i?i), Li 2i3 (i?i,i? 2 ), 

(4.1) oLi 4 (fli) logi? 2 , U 2 , 2 {Ri,R2) logi? 3 , 

(3.2) oLi 3 (i?i)Li 2 (i? 2 ), 

(3, 1, 1) o Li 3 (i?i) logi? 2 logi? 3 , (5-36) 
(2,2,1) oLi 2 (i?i)Li 2 (i? 2 ) log #3, 
(2, 1, 1, 1) Li 2 (i?i) logi? 2 log #3 logi? 4 , 

(1,1,1,1,1) Ologi?i logi? 2 log ^3 logi?4 logi? 5 , 

where the Ri correspond to the rational functions as indicated in section 5.2. 

Since the proposition implies TixTx = 0, we get Hx(S — Sx) = 0. Therefore we 
are essentially working at each step in the quotient space J-^ / Fx an d have to solve for 
considerably fewer coefficients than if we worked in F^ . If at any stage we cannot find 
a T A + with the desired property, i.e. cannot solve for the corresponding coefficients using 
our limited spanning set of input functions , the algorithm stops (we can of course try 
to rerun it with a larger set of input functions). Otherwise, it ends with producing an 
integrable tensor S'(i v .. i i) with S — -S^i,...,!) = 0, solving the main part of our integration 
problem. 

5.6 Elements in the kernel of the symbol map 

The algorithmic approach we described in the previous section often allows us, given a 
tensor S, to construct a function F such that S{F) = S. Let us now assume that the 
tensor S was obtained in some way from an analytic expression Fq (representing, say, a 
Feynman integral). It would be premature to conclude that the function F we constructed 
is equal to the original expression Fq, because they are only equal up to terms that are in 
the kernel of S. 

In the following we describe a way that, at least in most of the cases we studied so 
far, allows to fix this remaining ambiguity by parametrizing F — Fq in a suitable way. The 
parametrization we are proposing takes the form 

F-F = 5>^ + £ £ c h .., t k h .., e ■ ■ ■ , (5.37) 

W\ + ...+W(,=W 

where b\™ k ^ are defined in the previous section, and Ci l ...i l and Cj are rational coefficients, 
and ki 1 ...i e and k% are generators of the kernel of the symbol map. Below we give a (non- 
exhaustive) list of such generators, which cover a wide range of applications. The free 
coefficients can then be fixed by considering special values for the variables Xi, e.g., xi = 
or Xi = 1, yielding a linear system for the coefficients. 

In order for the above procedure to work we need to know the generators of the kernel 
of S. Even though this is a very difficult question to answer in general, we can compile a 
list of (presumably transcendental) numbers whose symbol should be defined as zero. 
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• <S(log(— 1)) = 0. An associated polygon would be "^J^' , and this gives the symbol 

o 

— 1 which is zero (modulo 2-torsion which we ignore). In particular, this shows that 
the symbol does not detect the multivaluedness of the logarithm (let alone of the 
polylogarithms), and in particular does not suggest how to fix the branch cuts of the 
polylogarithms. 

• All multiple zeta values (MZV's) are in the kernel of S. Indeed, MZV's can be defined 
as the values in X{ = 1 of the multiple polylogarithms Li mi) ... jmfe (xi, . . . , x&). Then it 
is easy to see that the associated decorated polygon will have all decorations equal 
to or 1, and hence the symbol vanishes. 

• In addition, there are combinations of transcendental numbers that individually have 
a non- vanishing symbol 5 , but there is a linear combination with zero symbol, e.g., 

5 ( Li4 (^) + 2l l0g42 ) = "( 1_ ^^^ + 2l (2020202) = O - (5 ' 38) 

The previous example is a special case of a more general result for so-called colored multiple 
zeta values, defined by the alternating sums 

s ni s n2 ■ ■ ■ s nk 

((m 1 ,...,m k ;s 1 ,...,s k )= ^'.^ ZmT > ( 5 - 39 ) 

I Li I ley ' ' ' I At, 

0<ni<n 2 <---<n fe K 

with rrii G N and G {±1}- It is easy to check that £(1, 1, 1, 1; —1, —1, 1, 1) = — Li4 (|) , 
and so eq. (5.38) can be written as 

S(((l, 1, 1, 1; -1, -1, 1, 1) - ^ log 4 2) = . (5.40) 

More generically, we have, 
Proposition 4. 

1. If at least one of the rrii is different from ±1 and (mi, s\) / (1, 1), then 

5(C(m,i, . . . ,m fe ; si, . . . , s k )) = . (5.41) 

2. ysi £ {-1,1}, Vm > 1, 
5 I CO^l; -1, s 2 ,...,s m )- ^ log m l - I -U. (5.12) 

m times 

The proof of this proposition will be given in appendix C. 



3 We recall that the "refined" symbol of a constant is not necessarily zero. 
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6. Application: a spanning set for harmonic polylogarithms 



In this section we illustrate our approach by expressing all harmonic polylogarithms (HPL's) 
up to weight four in terms of a spanning set of functions. We start by determining the 
arguments our spanning set should depend upon. Indeed, in the case of harmonic polylog- 
arithms we can classify all the arguments of the spanning set of functions under the mild 
assumption that all the arguments should be rational functions. Then, it is easy to see that 
the polygons associated to HPL's correspond to polygons where the root edge is decorated 
by the variable x, and all other sides are decorated only by or ±1. This implies that 
the tensor associated to an HPL is an element of the tensor algebra of the vector space 
generated by the formal basis vectors [x], [1 — x], [1 + x] and [2]. Indeed, the decorations of 
the polygon associated to an HPL are all ±1 or 0, except for the root which is decorated 
by the variable x. It is then easy to see that any dissection of this polygon will only involve 
the following five bigons 

n(m=^j=x, ii( ^==? J = 1 - x , 
o 1 

K^fe?) = 1+x ' K^) = K^M =2 - 
-i i -i 

Hence, the sets Phpl and Phpl of irreducible polynomials defined in section 5 are 

Phpl = {2, x, 1 - x, 1 + x} and P HPL = P HPL U {2 ± x, 3 ± x, 1 ± 2x} . (6.2) 

The most general rational function we can construct out of the irreducible polynomials in 
Phpl then reads 

±x ai (l-x)° 2 (1+x)" 3 2 a4 (2-x)° 5 (2+x)" 6 (3-x)" 7 (3+x)° 8 (l-2x)° 9 (l+2x) Ql ° , (6.3) 

with aj G Z. In the case of HPL's, however, it turns out (a posteriori) that we can restrict 
ourselves to the following set of rational functions 

1Z HPL = {±2 s x a (1 - xf (1 + x) 7 | a, 0,7,5 G Z} . (6.4) 

In the following we denote the elements of TZhpl by 

R% lS {x) = ±2 5 x a {l-xf{l + x)\ (6.5) 



with a, f3,j,5 G Z. It follows then from the previous section that we should consider 
the subset of these rational function: 
eq. (5.8)) that satisfy the constraint 



the subset of these rational functions (those contained in the set 7t HPL ; see eq. (5.7) and 



1 



for some integers a', (3', 7' and 5', and s = ±1. A little algebra shows that quadruples 
(a, /3, 7, 5) are confined to the values given in Tab. 2. The first line corresponds to constant 



- 30 - 



Q 


LZ 


R 
P 


7 


A 
c 






yJL 


P 


7 


u 


















-1 


+ 











-1 




1/2 


+ 


1 











X 




1 













—x 


+ 





1 








1 — X 


+ 








-1 







1/(1 + x) 


+ 


2 











2 

X 


+ 





1 


1 







1 - x" 1 




i 
l 


-1 


U 


U 


x/(x — 1) 


+ 


1 

1 


U 


-i 

-l 


U 




x/(x + 1) 


+ 





1 


-1 





(l-x)/(l + x) 







1 


-1 







(x-l)/(x + l) 




2 


-2 


-2 





x 2 /(x 2 - 1) 


+ 





2 


-2 







(1 -x) 2 /(l + x) 2 


+ 





1 





-1 


(l-x)/2 


+ 








1 


-1 




(l + x)/2 




1 


-1 





1 


2x/(x - 1) 


+ 


1 





-1 


1 




2x/(x + 1) 


+ 


2 





-2 


2 


4x/(l - x) 2 




2 


-2 





2 




-4x/(x + l) 2 



Table 2: Solutions to the constraint (6.6). Only half of the solutions are shown, all other solutions 
being related by i?f_ a)( _^ )( _ 7)( _ 5) (a;) = l/R s aM5 {x). 



arguments, and will not be discussed any further. Note that we could also include the 
inverses of the arguments in Tab. 2. Using the inversion formula for the classical polyloga- 
rithms, we can however always express Li n functions with inverted arguments in terms of 
polylogarithms taken at the arguments in Tab. 2, 

Lin ^— ^ = (— l)" _1 Li n (x) + products of lower weight terms, n > 2 . (6.7) 

Since furthermore the arguments in Tab. 2 are less than 1 for x G [0, 1], we will in the 
following only consider these functions as arguments of the Li n functions. Note however 
that even the functions in Tab. 2 are not completely independent, because we can use the 
distribution formula for the Li n function to relate three of them, e.g., 

Li n (x 2 ) = 2 n ~ l (Li n (x) + Li„(-x)) , (6.8) 

and three others using the same equation with x replaced by We finally arrive at the 
conclusion that, if we want to reduce all HPL's to some small (possibly minimal) set of 
(multiple) polylogarithms, the Li n sector of that set contains classical polylogarithms with 
16 different arguments. As we will see below, for lower weights we can find more relations 
among the spanning set of functions, reducing its size even further. In the following we 
will also have to deal with multiple polylogarithms of depth greater than one. From the 
previous section we know that the pairs of arguments (i?i,i?2) of these functions come 
from Tab. 2, subject to the additional constraint R\ — R2 € T^hpl- 

Finally, using the results from the previous section, as well as some elementary identi- 
ties among classical polylogarithms, we find the following spanning set of indecomposable 
functions up to weight four, 

• for weight one, 

#i 1} (x) = logx, Bf\x) = log(l-x), Bf\x) = log(l + x), ^ 4) (x) = log2, 

(6.9) 
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• for weight two, 



B^\x) = Li 2 (x), B^\x) = Li 2 (-x), B?>(x) = Li 2 



j(2) 



(3), 



1-X 



(6.10) 



• for weight three, 



£ 3 1J (x) = Li 3 (x), ^ 2) (x)=Li 3 ( 



?( 5 )^ - 



+ x 



B^(x) = Li 3 [ T 



tff(x) =Li 3 



Bf\x) = Li 3 (l-x), 

1 + X\ R (6), X T • 

, Z^ ; (x) = Li 3 



2x 



1 - x 



(6.11) 



• for weight four, 



B?\x 
Bf{x 

Bf\x 
Bf(x 



Li 4 (x), Bf\x) = Li 4 (-x), 
Li 4 (l-x), 2^ 4) (x) = Li 4 

x 



1 



= Li 4 
= Li 4 
= Li 4 



x — 1 
fl + x 
V 2 
1 -x 
1 + x 



1 + x 



i3^(x)=Li 4 
Hf(x) = Li 4 



^ 10) (x) 



Li 4 



i^ 13) (x) = Li 4 (1-x 2 



£i 12) (x) = Li 4 



x + 1 
1 — X 
2 

x — 1 

X + 1 

2x 



^ 14) (x) = Li 4l ; | 



x - 1 

2 



X 



£i 15) (x) = Li 4 



4x 



(x + 1) 5 



(6.12) 



(6.13) 



These functions are sufficient to express all harmonic polylogarithms up to weight three. 
Starting from weight four, we need to extend the set of functions by adjoining multiple poly- 
logarithms. We find that it is enough to add the following three supplementary functions 
in order to express all harmonic polylogarithms up to weight four, 

BT\ X ) = Li 2 , 2 (-M), Bf\x) = Li 2 , 2 Q, , BT\x) = Li 2 , 2 (1, ^) . 

(6.14) 

Note that, if the vector of singularities only takes values in the set {0, 1}, we can restrict 
ourselves to the smaller spanning set, 



for weight one: B^\x), B^'(x), 



• for weight two: B^\x), 



-32- 



• for weight three: B^\x), B^\x), 

• for weight four: B^\x), B^\x), B ( ^(x). 

Our choice for the spanning set is of course not unique, and we might have chosen a 

(i) 

different set of functions, related to the Bf functions via functional equations. Our choice 
was motivated by the fact that B^\x) is manifestly real for x G [0, 1]. Note that outside 
this interval, the branching structure of these functions can be more complicated. This 
issue is addressed in appendix E. 

6.1 Example 

Let us conclude this section by giving an example of how we can apply the procedure of 
section 5 to express a generic HPL of weight four in terms of the functions B^ . We discuss 
in detail the example of H(0, 0, 1, 1; x) = £2,2 0*0, ah other cases being similar. For the 
list of all results, we refer to appendix G. We start by deriving the tensor associated to 
H(0, 0, 1, 1; x). The polygon associated to H(0, 0, 1, 1; x) = G(0, 0, 1, l;x) is 



X 




There is only one relevant maximal dissection of this pentagon (all other dissections give 
rise to twogons with decorations and / or 1 which vanish by definition), 




From this dissection we can immediately read of the symbol of H(0, 0, 1, 1; x) as 

S{H(0, 0, 1, 1; x)) = (1 - x) <8> (1 - x) <g> x <g> x . (6.15) 

Note that in general the symbol of a harmonic polylogarithm H(a\, . . . , a w ; x), with cij G 
{0, 1} is simply given by (— l) k (a w — x) ... (a± — x), where k is the number of a^s 
equal to 1. Before turning to the question of how to express H(0, 0, 1, l;x) in terms of 
the spanning set, let us review how we could have obtained the symbol (6.15) using the 
recursive definition of the symbol (3.24). Note however, that in this case we cannot apply 
eq. (3.23) immediately, as the arguments of G(0, 0, 1, 1; x) = H(0, 0, 1, 1; x) are not generic 
and we would hence obtain divergences in the right-hand side of eq. (3.23). We therefore 
need to use a regularized version of the differential equation (3.23) [2], 

dff(0,0,l,l;x) = #(0,l,l;x)dlogx, (6.16) 

and so 

S(H(0, 0, 1, 1; x) = S(H(0, 1, 1; x)) ® x . (6.17) 
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We can iterate this procedure to compute the symbol of H(0, 1,1; x). The differential 
equation for H(0, 1, 1; x) is 



dff(0,l,l;x) = #(1, l;x)dlogx, (6.18) 

and so we get 

S(#(0,l,l;x)) = <S(tf(l,l;x))®x. (6.19) 

The symbol of H(l, l;x) = \ log 2 (l — x) is easy to obtain from eq. (3.7), 

S(H(l,l;x)) = ^S(log 2 (l-x)) = -<S(log(l-x))III<S(log(l-x)) = (l-x)<g>(l-x) . (6.20) 

Putting everything together, we immediately arrive at the symbol given in eq. (6.15). 

Let us now turn to the actual problem we want to study, namely how to express 
H(0, 0, 1, 1; x) in terms of the spanning set for HPL's defined in the previous section. The 

(k) 

goal is to find rational numbers c) ,■ such that the tensor associated to 



18 8 4 3 

Y/P b? w + EE i 2) B f (x) b? (*) + E i 3) *P (*) *?> ^ 

i=l i=l j=l i,j"=l 

+ EE C S ^ (x) B?> (x) Bf (x) + £ (x) *?> (x) Bf> (x) B« (x) , 

i=l j',fc=l i,j,k,l=l 

(6.21) 



equals the tensor given in eq. (6.15). The symbol of eq. (6.21) is easily obtained from the 
fact that S is linear and maps products of polylogarithms to shuffles, i.e., 



18 8 4 

5>S 1} 5(B w (x)) + £ y: i 2) 5(4^)) ni $(*?>(*)) 

i=l i=l j=l 

+ Ei 3) 5(6f(^))ffl5(6?w) 

+ EE cg<S(^ ) (x))UI<S(^ ) (x))III<S(Hf ) (x)) 

i=l j,k=l 

+ E ^5(Hf ) (x))m5(^\x))m5(Bi fe) (x))m5(H?(x)) 

i,j,k,l=l 



(6.22) 
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(i) 

The symbol of each function B) can be easily obtained, e.g., 



X \ \ I X \ X \ IX 



X + 1 J J \ X + 1 J \x + 1 / \x + l J \ x + 1 
X \ ( x \ fx 

<8 



x + lj \x + lj \x + lj \x + l / 

= (x + l)<8>x<8>x<8>x — (x + 1) (g) x <8> x <8> (x + 1) — (x + 1) <8> x <8> (x + 1) <8> x 
+ (x + 1) <8 X <8 (x + 1) <8 (x + 1) - (x + 1) <8 (x + 1) <8 X <8 X 
+ (x + 1) <8 (x + 1) <8 x (8) (x + 1) + (x + 1) <8 (x + 1) <g> (x + 1) <g> x 
- (x + 1) (8 (x + 1) (8) (x + 1) (8) (x + 1) . 

(6.23) 

The different shuffles in eq. (6.22) can be distinguished further by acting with the projectors 
defined in the previous section. In particular, we obtain 

18 

n 4 5(F(0,0,I,I;x)) = ^c J (1) n 4 5(Bj ) (x)). (6.24) 

i=i 

Equating the coefficients of the different elementary tensors on both sides of this equation, 
we obtain a linear system that allows us to solve for the for the coefficients 1 < i < 18. 
The solution is easily obtained, and is given by 

c« = _<£> = C W = i , (6.25) 

and all other coefficients are zero. We now proceed recursively, and subtract the (inte- 
grable) tensor arising from the solution we have found. By construction, the symbol of this 
difference must vanish under the action of the projector II4, 

n 4 <S(#(0,0, l,l;x)-L(x)) = 0, (6.26) 

where we defined 

L(x) = Li 4 (x) - Li 4 (l - x) + Li 4 (zTT[) " ( 6 - 27 ) 

(2) 

We next turn to the determination of the coefficients c^ - , i.e., coefficients of terms of the 
formB { j\x)B[ j \x). We can isolate these terms by applying the projector 113(8)111, 



(n 3 <8 III) S (H(0, 0, 1, 1; x) - L(x)) 

= £ i 2) (n 3 ^ni)[5(^ ) (x))ui5(^ ) (x))] 

ifei (6-28) 
J 2 ) rrr «r(')u\i <> rrr cz-rCj)/ 



Ei 2) [n 3 5(^ ) (x))]®[n 1 5(^ ) (x))], 



and equating the coefficients of the elementary tensors on both sides of the equation we 

c[j\ We find c^f = —1, and cf^ 



(2) (2) (2) 

can solve for the coefficients c\, . We find c 12 = —1, and q- = if 7^ (1,2). We 
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again subtract this contribution to find an expression that vanishes under the actions of 
both II4 and II3 IIi. Next we act with the projector II2 II2 on this difference 



(n 2 U 2 )S (H(0, 0, 1, 1; x) - L(x) + Li 3 (x) log(l - x)) = , (6.29) 

(3) 

and we immediately conclude that c\- = 0, VI < i,j < 3 . Similarly, by acting with the 

projector n20ni0lli we conclude that all the coefficients cf^ k must vanish. The remaining 
terms are thus all associated to products of logarithms, which can immediately be read off 
from the tensor 

S (H(0, 0, 1, 1; x) - L(x) + U 3 {x) log(l - x)) 

= ^ (1 - x) (1 - x) (1 - x) (1 - x) - ^x (1 - x) (1 - x) (1 - x) (6 _ 30) 

= S ^ log 4 (l - x) - ^ \ogx log 3 (l - x^j . 

At this stage, we have found a combination of (a product of) functions in our spanning set 
that has the same symbol as H(0, 0, 1, l;x), and so the quantities are equal up to terms 
that are mapped to zero by S. We make an ansatz assuming we have found sufficiently 
many elements in the kernel of S, as 

ai (hu + ^ log 4 2) + a 2 vr 4 + C 3 J2 b * B i ] (*) 

+ \ Y J c i Bf{x)+Y J C ij B { l\x)Bt\x) \ , 
\i=i i,j=i J 

where ai, hi, Ci and Cij are rational numbers, and ("3 = C(3) denotes the value in s = 3 of 
the Riemann £ function 



00 ^ 



n 

n=l 

The coefficients can now be fixed by looking at particular values of x. In particular, 
harmonic poly logarithms are known analytically up to weight four for x = 0, x = ±1 and 
x = ±2- It turns out that in all cases these values are enough to fix all the free coefficients. 
In the case of H(0, 0, 1, 1; x), we finally arrive at 

ff(0,0,l,l;x) =S 2 ,2(aO 

= — Li4(l — x) + ~L\i{x) + L14 



x — 1 

I 'l (6-33) 

- Li 3 (x)log(l - x) + xrlog 4 (l - x) - -logxlog 3 (l - x) 

24 D 

1 7T 4 
+ C 3 log(l ~X) + -IT 2 log 2 (l + — . 
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7. Conclusion 

In this paper, we have provided a review of the symbol map, a linear map that associates 
to a multiple polylogarithm of weight n an n-fold tensor in a way that captures many of 
the combinatorial properties of poly logarithms, and also respects the functional equations 
they satisfy. While so far the symbol map was defined recursively via iterated differentials, 
we have given a diagrammatic rule where the symbol of a multiple polylogarithm is ob- 
tained directly via a weighted sum over the maximal dissections of the decorated polygon 
associated to the polylogarithm introduced in ref. [21]. 

Furthermore, we have addressed the problem of integrating a symbol, i.e., the problem 
of finding a function whose symbol matches a given tensor that satisfies the integrability 
condition (3.15). We have presented a systematic approach of how to find a candidate 
spanning set for such a function. Once this candidate spanning set has been constructed, 
and working under the assumptions that its elements suffice to express the integrated 
symbol in, we showed how a set of projectors can be defined which help to find a function 
whose symbol matches the initial tensor. While our approach falls short of a complete 
algorithmic proof and is surely not adequate in all possible scenarios, we nevertheless 
believe that it can be applied in many situations, as was for example demonstrated in 
ref. [71, 72, 73] where our method was successfully applied to obtain new compact analytic 
results for certain one-loop hexagon integrals in D = 6 dimensions. Finally, we have used 
our approach to derive a spanning set for harmonic polylogarithms up to weight 4, and this 
spanning set was recently used to obtain an efficient numerical implementation of harmonic 
polylogarithms up to weight four [74]. 
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A. Review on shuffle algebras 

As shuffle algebras are a recurrent theme when working with multiple polylogarithms and 
their symbols, we review in this appendix the most important notions. Before turning to 
the special case of shuffle algebras, we first review some basic notions about algebras in 
general. 

Algebras over a field. An algebra A over a field F (F = R or F = C, say) is a vector 
space 6 over F together with an associative and distributive multiplication A <8> A — > A. In 

6 More generally, we could consider A to be a module over a ring R. 



-37- 



the case the multiplication has a unit element, the algebra is called unital. Furthermore, 
an algebra is said to be graded if A can be written as a direct sum as a vector space, 



A = @A n , (A.l) 

nel 

and if Va G A m and V6 G A n , we have 

a • b G An+n . (A.2) 

One of the most prominent representatives of a graded algebra is the tensor algebra asso- 
ciated to an F- vector space V, defined by 

oo 

T(V) = QT n {V), (A.3) 

n=0 

where %(V) = F and T\{V) = V, and for n > 2 we define 

T n (V) = V®...®V. (A.4) 

n times 

The multiplication on T(V) is defined on elementary tensors a\ ® . . . <%> a n by 



(ai (8> • • • <8 a m ) <8 (61 <8 . . . <8 b n ) >->■ ai (8) . . . <8 a m (8) b\ (8) . . . <8> 6 n , 



(A.5) 



making the tensor algebra into a graded algebra in an obvious way. Furthermore, the tensor 
algebra is also unital, the unit being the unit 1 G F. 

A homomorphism between two algebras A and B is a linear map that preserves the 
multiplication, i.e., a linear map <j> such that Va, 6 6 -4, ^>(a • b) = 0(a) • 4>{b). 

Ideals in algebras. An ideal in an algebra A (or more generally in a ring) is an additive 
subgroup X of A such that 

VaeA,Vbel, a-bel and b-a£X. (A.6) 

An easy example of an ideal is given by considering the ring of integer numbers Z and 
its subset raZ, n G Z, i.e., the set of all integer multiples of n. The set nL is obviously 
an additive subgroup of Z, and every time we multiply an element of nL by an integer 
number, we obtain another multiple of n. Hence nL is an ideal of L. Another example of 
an ideal is the kernel of an algebra homomorphism <f>. Indeed, Ker</> is a sub- vector space 
of A, and hence an additive subgroup. Furthermore, Va G A and Mb G Ker </>, we have 

cj)(a -b) = cj){a) -<f>(b) = 0(a) -0 = 0, (A.7) 

and so a ■ b G Ker 0, making Ker <f> into an ideal in A. 
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Shuffle algebras. After this rather general discussion on algebras, let us from now on 
focus exclusively on the example of shuffle algebras. As a starting point, let us consider 
a set £, whose elements we will refer to as letters, and consider the set W of all words 
constructed from elements in C, i.e., the set of all possible concatenations of letters in £, 
together with the empty word e, consisting of no letters (more precisely, W is the free 
monoid generated by the elements in C). We can define a multiplication on W given by 
concatenation of words, the empty word being the unit element. 

Let us now consider the vector space A over some field F given by all formal linear 
combinations of words in W. A is then in fact an algebra over F, the multiplication given 
simply by the concatenation of words. Furthermore, it is easy to see that A is also graded 
by the length of the words (the concatenation of two words with length m and n gives a 
word of length m + n). 

We can define another multiplication on A, the so-called shuffle product, defined on 
words by 

(ai . . . a m ) III (a m+ i . . . a m+n ) = ^ Oo-i(i) • • • a o- 1 (ni+n 2 )' (A- 8 ) 

ireE(ni,n2) 

where E(ni,n2) denotes the set of all shuffles of n\ + ri2 elements, i.e., the subset of the 
symmetric group S ni+ri2 defined by 

S(ni,n 2 ) = {a G S ni+Jl2 \ cr _1 (l) < . . . < cr _1 (ni) and <7 _1 (ni + l) < . . . < a~ 1 (ni+n2)} . 

(A.9) 

The vector space A together with the shuffle product is called a shuffle algebra. A shuffle 
algebra is again graded by the length of the words, and the empty word e is the unit 
element of the shuffle algebra. The definition (A. 8) of the shuffle product is equivalent to 
the following recursive definition, Vx, y G C, Vu, v G W, 

e III u = uYH e = u , 

(A.10) 

(xu) III (yv) = x(u III (yv)) + y{{xu) III v) . 

Note that the tensor algebra of a vector space V can be equipped with a shuffle product 
in a natural way: the set of letters C is simply a basis of V, the set of words corresponds 
to the elementary tensors a\ ® . . . <8> a n , the empty word is simply the scalar 1, and the 
concatenation of words corresponds to the multiplication (A. 5) of two tensors. 

In section 2 we have encountered another example of a shuffle algebra, the shuffle 
algebra of multiple polylogarithms, eq. (2.3). In that case letters are the elements 
of the vector of singularities (a±, . . . ,a n ), the latter being the words. Concatenation is 
simply defined by the concatenation of the vectors of singularities. The length of a word 
corresponds to the weight of the polylogarithm, i.e., the number of components of the 
vector of singularities. In other words, the shuffle algebra of polylogarithms is graded by 
the weight. 

B. Selected examples of symbols 

In this appendix we compile a list of the symbols of the most commonly used (multiple) 
polylogarithms. First of all, the symbol of an ordinary logarithm is simply the argument 
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of the logarithm, 

S(\ogx)=x. (B.l) 
From eq. (3.7) it follows then, for every non-negative integer n, 

S (Jj log" x] = x ® . . . ® x = x® n . (B.2) 

n times 

The symbol of the classical polylogarithms have a similarly simple form, i.e., for every 
non-negative integer n we obtain 

S(U n {x)) = -(1 - a;) <g> x® (n-1) = -(1-x) (8> a® (gig , (B.3) 

(rt— 1) times 

where we used the notation of eq. (B.2). Note that for n = 1, the classical polylogarithm 
can be expressed as an ordinary logarithm, Lii(x) = — log(l — x), which is consistent with 
the symbols given in Eqs. (B.l) and (B.3), 

S (Lix (X)) = -(1 - x) = S(- log(l - x)) . (B.4) 

Finally, the symbol of a Nielson polylogarithm reads, 

S(S n , p (x)) = (-If (1 - x)® p <g> x® n = (-If (1 - x) <g> . . . <g> (1 - x) ® x ® . . . ® x . (B.5) 

S v ' V * ' 

p times n times 

Again we note that the Nielson polylogarithms contain the classical polylogarithms as a 
special case, S n -i t ±(x) = Li n (x), an identity which is easily verified at the level of the 
symbols (B.3) and (B.5). 

The previous examples of the classical and Nielson polylogarithms are both just special 
cases of harmonic polylogarithms where the components of the vector of singularities take 
values in {0, 1}, 

Li n (x) = il(0 n _i, 1; x) and S UiP (x) = H(0 n , l p ; x) . (B.6) 

The symbol of a harmonic polylogarithm H(a±, . . . , a n ; x), with £ {0, 1}, can be written 
in the compact form 

S(H(m, ...,a n ; x)) = (-l) k (a n - x) <g> . . . <g> (ai - x) , (B.7) 

where k is the number of components in the vector of singularities (ai, . . . , a n ) equal to 1. 
Indeed, the polygon P(a n , . . . , a±, x) associated to G(a±, . . . , a n ; x) = (— l) k H(ai, . . . , a n ; x) 
has the root side decorated by x, and all other sides decorated by or 1. It is easy to see 
that the only relevant maximal dissection of such a polygon is 



X 
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All other maximal dissection give rise to bigons decorated only by and / or 1, which give 
a zero contribution to the symbol. The non-vanishing dissection produces a term in the 
symbol given by 



/i 



O-n-l 



(B.8) 



12 



ai 



which is equal to the tensor in the right-hand side of eq. (B.7) (apart from the sign). 

Generic harmonic polylogarithms where the components of the vector of singularities 
take values in { — 1,0, 1} do not admit a compact expression for the symbol. In this case, 
the symbol is however easily obtained from the symbols of generic multiple polylogarithms, 
which are reviewed up to weight four in the next subsections. 



B.l The symbol of a generic multiple poly logarithm of weight one 

A generic multiple polylogarithm of weight one can be written as G(a;x), with a, x G C. 
We can associate a bigon to this function via 



G(a;x) o P(a,x) 



(B.9) 



The symbol of G(a; x) is then 



S(G(a-,x))= fi (*=^) 



(B.10) 



B.2 The symbol of a generic multiple polylogarithm of weight two 

To a generic multiple polylogarithm G(a, b; x) of weight two we associate a trigon 



G(a, b; x) o P(b, a, x) 




(B.ll) 



The symbol of G(a, b; x) is then obtained by looking at all the maximal dissections of the 
trigon, obtained by inserting a single arrow. In the following we give the three maximal 
dissections, together with the term in the symbol they correspond to. We use the shorthand 



ab\cd = [x 



(B.12) 



The three maximal dissections of the trigon in eq. (B.ll), together with the term in the 
symbol S(G(a, b; x)) they correspond to, are 






+ax\ba 



+bx\ax 



-bx\ab 



B.3 The symbol of a generic multiple polylogarithm of weight three 

To a generic multiple polylogarithm G(a, b, c; x) of weight three we associate a tetragon 



G(a, b, c; x) o P(c, b, a, x) = 



(B.13) 



The symbol of G(a, b, c; x) is then obtained by looking at all the maximal dissections of the 
tetragon, obtained by inserting two non-intersecting arrows. In the following we give the 
twelve maximal dissections, together with the term in the symbol they correspond to. We 
use the shorthand 



ab\cd\ef = fi 



/i 



/ 



(B.14) 



as well as the notation for shuffles, 

A\BU1C = A\B\C + A\C\B . 



(B.15) 



The twelve maximal dissections of the tetragon in eq. (B.13), together with the term in 
the symbol S(G(a, b, c; x)) they correspond to, are 



b b b b 

+cx\bx\ax +ax\ca\ba — bx\abUlcb +cx\ac\bc 





b b b b 

+cx\ax\ba +ax\ba\cb +cx\bc\ab — cx\axYUbc 



b b b b 

— cx\bx\ab — ax\ca\bc bx\axUlcb — cx\ac\ba 

B.4 The symbol of a generic multiple polylogarithm of weight four 

To a generic multiple polylogarithm G(a, 6, c, d; x) of weight two we associate a pentagon 



G(a, b, c, d; x) -H- P(d, c, b, a, x) 



_ d 




(B.16) 



The symbol of G(a, b, c, d; x) is then obtained by looking at all the maximal dissections of 
the pentagon, obtained by inserting three non-intersecting arrows. In the following we give 
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the 55 maximal dissections, together with the term in the symbol they correspond to. We 
use the shorthand 

ab\cd\ef\gh = n(»=^J ® //(«^L^ ® M («=L^ ® j (B.17) 



as well as the notation for shuffles 



A\B\CUID =A\B\C\D + A\B\D\C , 
A\BUI(C\D) = A\B\C\D + A\C\B\D + A\C\D\B , 
A\BU1CU1D = A\a(B)\a(C)\<T(D) , 



(B.18) 



where in the last equation the sum runs over all permutations of the set {B, C, D}. The 55 
maximal dissections of the pentagon in eq. (B.16), together with the term in the symbol 
S(G(a, b, c, d; x)) they correspond to, are 








+dx\cx\bx\ax +ax\da\ca\ba — bx\abUl (db\cb) +cx\dcUI (ac\bc) — dx\ad\bd\cd 




+cx\dcUI (bx\ax) — dx\ad\ca\ba +bx\ax III (db\cb) —cx\dcUI (ac\ba) +dx\ad\bd\cb 




+ax\da\bd\cd — dx\cx\bx\ab — ax\da\ca\bc +bx\abUl (db\cd) +dx\cx\ac\bc 




— cx\dcUI (bx\ab) +dx\ad\ca\bc — bx\ax III (db\cd) — dx\cx\ac\ba — ax\da\bd\cb 
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+bx\ax III (cb | dc) +dx\cd\ac\ba — dx\axUI (bd\cb) +cx\dcUI (ax\ba) — dx\ad\ba\cb 




— ax\ba\db\cd +dx\cx\bc\ab +ax\da\cd\bc +dx\cd\H(bx\ab) — dx\ax\ca\bc 




+ax\ca\dcUI ba +dx\bd\abUl cb — cx\dcUl 6c III ax +dx\ad\cdUl ba +dx\bx\ax III cb 




—ax\ca\dcUIbc — dx\bd\cdJJI ab — dx\cx\ax III be — ax\da\cdJJI ba —dx\bx\cbUlab 




+ax\ba\cb\dc — dx\cd\bc\ab +dx\ax III (cd\bc) — dx\cd III (ax\ba) +dx\ax\ba\cb 
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C. Proof of Proposition 4 



In this section we present the proof of Proposition 4. The proof uses the combinatorics of 
the decorated polygons introduced in section 3. In order to be able to map the symbol of 
colored multiple zeta values (CMZV's) to polygons, we first have to relate the CMZV's to 
multiple polylogarithms. From the series representations (2.9) and (5.39), it is easy to see 
that one has the relation (provided that the CMZV's converge) 

((m 1 ,...,m k ;si,...,s k ) = {-l) w G mi ,..., mfc (li, . . • , s k ) , (C- 1 ) 

where we defined w = mi + . . . + m k and 

3 



Us*. (C.2) 



i=l 

Hence, using the correspondence between multiple polylogarithms and decorated polygons, 
we can associate to the CMZV ((mi, . . . , m k ; si, . . . , s k ) the polygon 7 

P( 0,...,0 0^0 ,si,l). (C.3) 

m^ — l times mi — 1 times 

We start by introducing concepts needed to prove, and then prove Proposition 4. We 
will in fact proof results at the level of the polygons that are slightly more general than 
the results given in Proposition 4 but not imposing the restriction (mi, si) 7^ (1, 1) (which 
corresponds to divergent CMZV's). The first proposition we give is a generalization of 
statement 2 in the proposition. Statement 1 in Proposition 4 is equivalent to Proposi- 
tion C.6. 

Proposition C.l. The symbol of P(ei, s n , 1) for some e; L = ±1 is equal to \ a ,n(2® n ) 
for 

A a ,n = (-l) a ^ n ^ and a = n - max{i | Si = -1} . (C.4) 

We start by noting that from the Holder convolution (2.8) with p = 00 it follows that 
P(xi, x m , 1) has the same symbol as the polygon P(l — x m , 1 — xi, 1). So, without 
loss of generality, we consider the polygon P(0, 0, 2, 0, 0, 2, 0, 0, 2, 0, 0, 1) and 

to ti t m 

will find its symbol, remembering to take into account the factor of (— l) n . The move 

from sides labelled 1 and — 1 to sides labelled and 2 increases the number of dissections 

that have coefficient 0. The combinatorics of the dissections of polygons of this type is 

best captured by a certain type of planar trees, the so-called Hook-arrow trees, which is a 

change of view on the maximal dissections, and hence the symbol, of a polygon. 



7 There is of course the factor (—1)™ of eq. (C.l) to be kept in mind. 
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Hook-arrow trees. Every full dissection of a polygon uniquely defines a certain spanning 
tree r on the vertices which are the midpoints of the polygon sides, and vice-versa. These 
vertices, v\,...,v n , inherit the label of the side they sit on and they form the vertices of 
r. We induce the edges of r as all possible lines, between the Vi, that do not cross arrows 
from the dissection. Here is an example of a full dissection of a 4-gon with the spanning 
tree induced: 



4 




2 



We also induce a root on r as the vertex lying on the final side of the polygon. The 
edges are then oriented towards the root. For the above example of a dissected 4-gon the 
rooted spanning tree is: 

4 




2 



The edges of the tree will not cross by construction; we define interlacing to reflect 
this for use in the definition of a hook-arrow tree. 

Definition C.2. A tree with a linear order on its vertices Wj is said to be interlaced 

if there exists a choice of four vertices w± < ... < W4 such that both edges {wi,ws} and 
{^2,^4} are contained in the tree. 

We now give a formal definition of a hook-arrow tree and illustrate the definition with 
an example. 

Definition C.3. A hook-arrow tree is a rooted spanning tree on a set of vertices in a 
linear order, v\ < ... < v n , which is not interlaced and has root v n . 

We can think of a hook-arrow tree as being a tree embedded in the plane on the vertices 
arranged in a circle. 
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Example. For the polygon P(2, 0,2, 0,0, 2, 0,1) (attached to the multiple polylogarithm 
G (0,2, 0,0, 2,0,2;!)) we have the following possible maximal dissection 



1 




We note that, for clarity, in this example we give the arrows in the dissection of a polygon 
dashed shafts. 

We now construct the hook-arrow tree for this dissection, following the method outlined 
in Definition C.3 (see also fig. 1): 

1. Add vertices to the middle of each side of the polygon; 

2. Join all vertices that can be connected without crossing the shaft of an arrow; 

3. Remove the polygon and arrows and direct tree towards distinguished vertex repre- 
senting final side of polygon. 



1. 2. 3. 

1 1 1 




Figure 1: Construction of a hook arrow tree corresponding to a maximal dissection. Details are 
given in the text. 



We now reintroduce the dual tree view of a dissection from section 3.1 of ref. [21] as 
this is also beneficial in finding the symbol attached to a polygon. As with the dissection 
of a polygon using arrows, the dual tree can easily be seen in the hook-arrow tree view. 
For clarity we give the dual tree a dash-dotted line. 
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1 1 




Hook-arrow tree and dual tree Polygon dissection and dual tree 
We now elaborate on the interpretation of a hook-arrow tree. 

The association between a 2-gon in a polygon dissection, an edge of a hook-arrow tree 
and a term of a tensor product in the symbol is 

Term represented 



i — 

a 




In the dissection of a polygon, we set the coefficient of certain terms in the symbol to 



if the sides of 2-gons have certain combinations of labels, e.g., [i 



0. We do 



the same in a hook-arrow tree when we have the corresponding edges 







for any a. 

To obtain the term in the symbol, the edges of the hook-arrow tree are chosen in the 
same order as the corresponding 2— gons in the dissection of a polygon. This can also be 
seen by viewing the dual tree (this correspondence is shown in the above example). Note 
that an explicit algorithm for this purely from the view of the hook-arrow tree can be 
written. 

The sign of a maximal dissection is determined, in a similar way to the polygons by first 
determining the number, a, of 'backward' edges in the hook-arrow tree. These correspond 
to edges that, respecting their direction, go from one vertex in the order of the vertices 
defined by the polygon, to a previous vertex. The sign is then (—1)". 

Remark C.4. As previously mentioned, it is important to note that for every full dis- 
section of a polygon there is a unique corresponding hook-arrow tree and that the method 
of extracting the symbol from hook-arrow trees is simply a different view of the procedure 
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for polygons. We include an overview of the construction of an hook-arrow tree because 
constructing all possible dissections for the functions required in Proposition C.l is a lot 
easier to view from this perspective. The actual proof then takes place in the dual tree view. 



Finally, before the proof of Proposition C.l we require the following proposition for 
which we give a sketch proof involving generating functions. 

Proposition C.5. If c, n G Z-° then 

g ( - 1)l Cn-t C )( n + ' +1 ) =( " ir 

Proof. (Sketch) 

Let r = n—i and view the right hand side of the identity as coefficients of the generating 
function 

*(*)=£*»B-ir-( r+c )(" +c+1 ). 

iA 7i \ r )\ n-r ) 

Firstly since for r > n we have ( n ^^ 1 ) = 0, we can change the summation of r to run over 
all positive integers. Then by re-ordering the summation signs for small x, and applying 
some basic properties of binomial coefficients we can arrive at 



Then by using 

m. x ' v ' m. x ' 



n+C+1 



x m = 



(1-X) 



k+1 



we can show that 

oo 
n=0 



n=0 

which proves the result. □ 



Proof, (of Proposition C.l) After applying the Holder convolution with p = oo, and without 
loss of generality, we attempt to find all hook-arrow trees relating to the polygon 

P((W), 2,(W), 2,0, ... ,0,2,(W),1) 

to tl t m 

which do not represent terms with coefficient in the symbol. After some consideration 
we see that these hook-arrow trees must take the following form. 
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Each ti t i and U^, for i = l,...,m — 1 are chosen integers < ii,i,ii,2 < U such that 
ti,i + tj,2 = U. The choice of the ijj arises from the fact that we can choose where to 
partition each group of tj vertices labelled 0, for i = l,...,m — 1, and attach them to 
the vertices labelled 2, remembering that the vertices must not cross. In the case of the 
function P(2, 0, 2, 0, 0, 2, 0, 1) from the example above, where m = 3, to = 0, t\ = 1, £2 = 2 
and <3 = 1, we would have 6 possible valid dissections, arising from two choices of the tjj 
in + ii 5 2 = 1 and three choices from £2,1 + *2,2 = 2. We note that the example above 
explored the particular dissection where t\ t \ = 0,ii ; 2 = 1,^2,1 = 1 and £2,2 = 1- 

We will now show how it is possible to simplify this tree by effectively removing the 
edges joining vertices labelled 2 and 1 and replacing them with edges connecting vertices 
labelled and 2. For this we return to the dual tree notation described above, and in 
section 3.1. The dual tree of the above hook-arrow tree above is 
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We will now write the tensor product of the symbol of the left hand dual tree part in 
the above claim. Let us recall out convention (3.9) that shuffles takes precedence over a 
tensor 

a® 6 = a(g> ...(g) a and that a® 6 III c = (a <g> ... <g> a) III c. 

V v ' V v ' 

b b 
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The left hand dual tree part in the claim will have the symbol 
(-l)'* 1 ' 1 Q <g> 2®**" 1 - 2 III III Q <g> 2® ifc - 2 III 2® c ^ 



*fc,i=° 

tm— 1 



*m-l,l=0 



^m— 1, 



r (_l)*fe+i+ 1 ^1 (g) 2®**-i.2 TTT 2®( c+tfc+1 )^ 



which is exactly the symbol for the tree on the right side. Note that we used Proposition C.5 
in the last line of the calculation. 

By repeated application of this simplification, starting with c = t m and k = m — 1, we 

m 

will arrive at a much simplified tree. By noting that n — a — 1 = ^^{U + 1) and recalling 

i=i 

that to = a we see that this tree is 




times a factor of (— l)( n a l \ This represents the symbol 



n — 1 



Finally, by including the factor of ( — l) n from the application of Holder involution we find 

□ 



Proposition C.6. The polygon P(0, 0, £i, . . . , 0, 0, 1) /or £j = ±1 and at least one 



mi— 1 



m k -l 



of the rrii / 1, /ias a symbol equal to 0. 
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Proof. (Sketch) After applying Holder involution (2.8) with p = oo we try to find possible 
hook-arrow trees which do not correspond to terms with coefficient in the symbol. The 
vertices of the hook-arrow tree will be labelled corresponding to the sides of the polygon 

P(li,i, -, 7to,i) 2 > 7i,2, -,7ti,2, 2, ...,2,7i jm , ...,7t m , m , 1). 

where all the 7^ are equal to either or 1. As in the proof of Proposition C.l, the vertices 
labelled 2 must connect directly to the final 1 and the vertices labelled must connect to 
a 2. However, there is no way to connect the l's to any other vertex without setting the 
coefficient of the term to 0. There is therefore no possible non-trivial dissection, and hence 
the symbol is zero. □ 



D. Some considerations on the implementation of the algorithm 

The algorithmic approach presented in section 5 relies on the construction of the sets TZ^ 
defined in eq. (5.17). The algorithm then consists in selecting the elements of TZ which 
satisfy certain factorization properties. Even though this is a mathematically well-defined 
prescription, implementing this algorithm into a computer program can be hampered by 
several issues, 

1. The set TZ is infinite, and so we cannot just proceed by 'trial and error' to select the 
elements that have the right factorization properties. 

2. Factorization of polynomials is rather slow (at least on most computer algebra sys- 
tems), leading to serious speed issues. 

The first issue can be dealt with by decomposing TZ as 

00 

K=\JlZ n , (D.l) 

n=0 

where 1Z n is defined as the subset of 1Z consisting of those elements ± nf=i sucn that 
\n\\ + . . . + \nk\ = n. In order to construct the set TZ^\ which is the basis out of which 
TiP*) for arbitrary k is constructed, one can then limit oneself to truncating the tower of 
sets in eq. (D.l) to some finite value N < 00. Indeed, in practical applications one does 
not expect rational functions for large values of the sum n of the exponents 8 . 

The second item seems to be harder to solve, since it is related to the capabilities of 
the chosen computer algebra system. It is however possible to circumvent this problem by 
deriving from the factorization constraints a necessary condition that must be fulfilled and 
that can be checked in a fast way by a computer. Since in practice most of the elements 
of TZ will fail this constraint, one can filter out these elements and discard them in an 
easier way. In the following we discuss the example of TZ^ l \ The generalization to TZ^ is 
immediate. 

8 Empirically, we observe that the set seems to be finite in general. 
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Let us consider a generic element R in 1Z. Without loss of generality, we can assume 
that we can write 



m — ti£ 

7Ti . . . 7T 



R = s l\ + r (D.2) 

Kl+l ■■■ 7r fc 

with s = ±1 and n« G N, and 7fj G P. Checking if i? G 72.W is equivalent to checking 
whether 1 — R G 72., i.e., whether 1 — i? can be written, up to a sign, as a ratio of elements 
from the set P. Writing 

— n( + i — riu — ni — no 
7T»li . . .IT, — STTi . . .IT/ 

1-R= + =4l di -> ( D - 3 ) 

^+1 ••• 7r fc 

it is easy to see that this condition can only be fulfilled if the polynomial in the numerator 
can be factored into a product of elements in P. Let us call LT this numerator. A necessary 
condition for R G TZ^ is thus that II can be divided by at least one element in P. We can 
further simplify this condition by reducing it from a problem of division of polynomials to 
a problem of division of integers. Indeed, we can choose prime numbers {pi, . . . ,p m } such 
that 7fj(pi, . . . ,p m ) ^ Tfj(pi, . . . ,p m ), for i 7^ j. The necessary condition for R G TZ^ then 
reduces to checking that the integer number H(pi, . ■ ■ ,p m ) can be divided by at least one of 
the numbers 7fi(pi, . . . ,p m ), which is in general much quicker to test on a computer. Note 
however that this is a necessary, but not necessarily sufficient, condition for R G TZ^\ 



E. Analytic continuation of the spanning set of functions for harmonic 
polylogarithms 

E.l Analytic representation inside the unit disc 

The analytic expressions for the spanning set {B^\x)} introduced in section 6 are valid 
for x G [0,1], but the functions might have different analytic representations in different 
regions of the complex plane. In particular, since the Bf functions are real when x is in 
the range [0, 1], Schwarz's reflection principle implies that the elements of the spanning set 
must satisfy 

Bf {x*) = Bf (x)*, (E.l) 

where x* denotes the complex conjugate of the complex number x. We checked numerically 

that the analytic expressions for the elements of the spanning set are valid everywhere inside 

(13) 

the unit disc, except for B\ (x), where the correct expression for \x\ < 1 is 

(13) _ j Li 4 (l - x 2 ) , if Re(x) > or (Re(x) = and Im(x) > 0) , 

4 {X) ~ \ Li 4 (l - x 2 ) - f a(x) log 3 (l - x 2 ) , otherwise , 

(E.2) 

where a(x) = sign(Im(x)). In order to understand this structure, let us look at the simpler 
case of weight two. Then we get 

7T 2 

Li 2 (l - x 2 ) = - log(x 2 ) log (1 - x 2 ) - Li 2 (x 2 ) + — . (E.3) 
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The first term in this expression contains log(x 2 ) = 21ogx, which is real for x > 0, but 
develops an imaginary part for x < 0, displaying the complicated branch cut structure of 
Li2(l — x 2 ). A similar reasoning leads to eq. (E.2). 



E.2 Analytic representation outside the unit disc: inversion relations 

We now have analytic representations of the elements of the spanning set everywhere inside 
the unit disc, and so we can analytically continue them outside the unit disc via inversion 
relations, i.e., functional equations of the form 

B^\x) = ^2 Cijki&k ( ~ } + products of lower weight. (E.4) 

k,l W 

These functional equations can easily be obtained for the whole spanning set. Below we 
show the explicit inversion formulas for weight one and two. For the the complete list of 
inversion formulas for higher weights, we refer to appendix F. Letting <j{x) = sign(Im(x)), 
we get, for \x\ > 1, x not real, 

• for weight one: 

ef'W =-<»(!). 

0? 1 to = -<> (i) + Sf (i) - ixc(x) , (E.5) 

S f'(^ S f. (I) -<»(!) , 

• for weight two: 

B <.) ( ^- W i"(i)-^(i) 2 -^(i)4, 



2 1 \x) \x) 6 

B f W - - >o g2S ™ (I) - 1 BS " (I) 2 + «P» (i) Bf' (i) + 3?' (1 

T 2 

4 ' 



(E.6) 



(2) / 1 \ . H (3) (l\ vr 2 



Note that the value of a is ambiguous for real values of x. This ambiguity can be resolved by 
the He' prescription commonly used in the physics literature. According to this prescription, 
we need to assign a small imaginary part to real value of x, i.e., if x is real, we need to 
perform the replacement x — > x + ie, and this replacement fixes at the same time the value 
of a. Note however that some care is needed when applying the inversion formulas to real 
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values of x because Schwarz' reflection principle implies 

B^\x ± ie) = ^2 c ijki — j^^^ P r °ducts of lower weight 

= ^2 c ijki B$ ( — T i£J + products of lower weight (E.7) 

k,l ^ X ' 

= ^2 c ijki B^ ( — ± ie] + products of lower weight. 

k,l \ X ' 



If x lies on the unit circle, \x\ = 1, there is an ambiguity whether to use the expression 
for Bj l \x) valid inside or outside the unit circle. We checked numerically that the two 

values agree in all cases, except for B^\x) = L14 (^ jj^yz ^J ■ In this case, we find 

i.e., the argument of L14 is real and greater than 1 for every x on the unit circle (except for 
x = — 1, where the result is divergent), and we have an ambiguity on the imaginary part 
of B^ (e l,p ) . This ambiguity can be lifted by requiring the function to be continuous in a 
neighborhood of the unit circle. To study this, let us consider a circle which is infinitesimally 
close to the unit circle, i.e., we choose x = (1 — e) e lip , for some infinitesimal e. We then 
find 

(TTx? = c^f (l + -tanf + (V)), (E.9) 
i.e., we see that for \x\ = 1, we have, 

,(15) , . - ■ ( 4x 



with a(x) = sign(Im(x)). 



B\ L0 > (x) = Li 4 ( j^—y + iue ) > (E.IO) 



F. Inversion formulas for the spanning set 

In this appendix we present the inversion formulas for the spanning set for weight three 
and four, valid for x £ C x . 

F.l Weight three 



^<*>-i4 ,, (;)' + i^ ,, (;) + <; 



(F.l) 
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^M^© -^'(i)^(i) -^(l)-* 



■(2) 



fif 1 (i) = imr(i) log2Bf ) (-) - im(x) log2Bf 1 ( - ] + 5 log 2 2E; 1 ' ( - 
-i 1 o^^(i) + i 1 o E <»(i) 2 - 1 o g2B f'(i) 8 !"(i) 



^-(S)/^ 3 1«(2) /l^ H (3) /^V , 1 j n (2)/l\ , 1 2„(3) 



- B'P (i) - Sf (I) + log 2 2 + i log 3 2 - i» 2 log 2 + C, 

fif (x) = ilog 2 2B« (I) - ilogM 2 ' (1) + ilo g2e (" (I)' 

^ w - 4'- 2 <' (i) - ^i 2 ' (;) + 1 ^! 2, (I) ^ (i) 

l-(3)/l\ 3 !„(2) K (3) /1\ 2 , 1 2«(3) „(5)/l 



4 6) fM+^og32-^ 2 lo g 2 + C 3, 
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fif W = i™(x) log 2tf > (i) + ilog 2 2fif» (i) - Il„g2B« 2 > (I)' 
- -»7r<7(s) log 2 2 - - log 3 2 + -tt 2 log 2 . 



F.2 Weight four 



^w-^»(i) 3 -^»(i) 4 + ^»(i) 2 



»4 W " 24 ^i [ x ) 12 * *i [ x ) *A y x ) 360 , 

(3) _ _i# V , Ir(2) A 1^ R (l) A 1 V _ Ir(2) / R (l) ^ 1 



8f (x)=Bf (i), 



45 ' 



tff (*) = B< 4 > 



fif) (x) = -^vra(x)log 2 2^ 1) (-) + ^tkt(x) log 2 2tff } f- 

- ii^( I )lo g 2S| 1 > (0 +! ™(x)log2Bf» B<" (i 

- l* roW ]og2B» (i) 2 -ilog 3 2B« (I) + ilo g >2Bf) (I 
-IV2 B «(i) 2 + IV2 B f)(I) B «(i)-JV2 8 «(I 



- 58 - 



1 .2, OR (3) ( 1^ 1 . , s H (l) 1 \ 3 , 1 ■ / n„(3) fl \ R (l) /" 1 



- H^" (i)V'(i) + >^ (Jf - 

+ (I) 2 _ e u, (I) _ i jlw(x) log3 2 _ >_ log4 2 + log2 2 + ^ 

fif W = - i log' 2S<» (1) + 1 log' 2S< 2 > (i) - \ log 2 2Bi" 

+ ^(;M;)'-^(;)>£ 

-Mi) 4+ ^(iMi)'-^i)>(i 
-^i) 4 -^<;)'-^(;H* a 

"12" bg 2 " 360' 



(x) = Bf> 



^ 10) (x) = B< 9) ( 1 



a; 

Bl 1 " (x) = - \>™(x) log 2 2B< 3> (i) + log2S? » (i V + i log 3 2S< 3) Q 

-i^ 2 ^G) 2 ^^ 2B ! 3, (i) 3 -k^ 2 ^ Vl 
-H^ 3, G) 3 -^! 3, (i) 4 + ^! 3, (i) 2 -^(i 

+ -ivra(x) log 3 2 - — log 4 2 + -vr 2 log 2 2 + ^ , 
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£i 12) (x) = ^ma{x) \og z 2B\ 



(2) 



log^ 2B\ 



(2) 



f -ma(x)B\ 



(2) 



-iira(x) log 2B\ 



(2) 



+ ^log 3 2^ 



+ -log2^ 2) 



1 B M 
24 1 



^ 2 log2Hf) 



B 



(8) 



i I I 7T 

-ivrcT(x) log 3 2 - — log 4 2 + -vr 2 log 2 2 + ^ , 



Bf > (x) 



"3^ 



(2) 



(i) 



+ 3^ 



(i) 



Bi 



(i) 



(2) 



2B 



(3) 



B[ 



(3) 



B\ 



(1) 



Bi 



(i) 



n 2 B[ l) 



Bi 



(i) 



+ Bi 



(2) 



+ Bi' 
1 

~ 24 



(2) 



Bi 



(3) 



Bi 



(2) 



Bi 

4 



(3) 



1 

X 

1 

X 

1 

24 



Bi 



(i) 



Bi 



(i) 



Bi 



(3) 



(1) 



+ ,- 2 B? ) 



Bi 



(i) 



Bi 



(i) 



B 



(14) 



Bi 



B\ 



(3) 



Bi 



(3) 



B 



(3) 
1 



Bi 



(3) 



7T 2 Bi 3) 



Bj 14) (x) = lma{x)B\ 



(2) 



+ -ma{x)B{ 



(3) 



77T 4 

360 ' 

2 



+ 2^(x)Bf } (1 ) 



(2) 



24 1 
+ ^ 2 Bf) 

24 1 



i) ^ 
- ^ 
- & 



l -ma{x)B { l ] 
6 



Bi 



(2) 



Bi 



(2) 



Bi 



(2) 



B 



1 

:r 

(13) 



+ W> 



I 



Bi 



(2) 



r 

X 



+ 



7T 



45 ' 
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(x) = Ama{x) log 2 2B^ ] Qj - Sma{x) log 2 2Bf ] Q 

+ 2i7r<r(:r) log 2^S 1} Q - 8i7ra(x) log 2fi[ 3) Q Q 
+ 8ma{x) log2^ 3) + ^trt^)^ (- 

-w W ^(i)^(i) 2 + w WB p(i) 2 B «(i 

-fi™(z)B< 3 > (i) 3 + B< 15, (i) + ^^(x)log 3 2, 



12 v ' 1 v^/ 2 W 24 1 v^/ 

-?^<>(i) 2 - s ™(i)^»(i) B p'(f 

- B ?'(i)^(i) , + ^ B ? , (i) > -s B ? , (i)> 

+ 5 log 3 2Bf (I) +B a»(i) B p»(i) 2 B (3»(I 

+ ^(i)v 2B f'Q)-i^(i) B p(i 

+ 2B f (i) B «(i) +B f-(I) 2 l„ g 2 B p(i 
+ ^log^Msf (I) - i»V(x) B f (i) + 3 C ,S?» (i~ 
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- 2S f(i)- 2B f(i)- 3S -(i) +3B f)(i)-? S f>(i 
- B r(i)- B r(i) + ^(i)io g2 -^»(i).o g ^ 

- ^log^fx) + i»> log 2«t(x) - jS! 1 ' (±) Ca + 6Li 4 Q) - ^ , 
Br » W = (I)' _ ^<.» (i) B p> (I)' _ ^ (i) B? (I)' 

-^o g ^'(i)^^ WB P(i) S -^(i) 2 ^(i) 2 
+ ^V 2B p(i) 2 + ^»(i) B P(i) B P(i) 2 



4 1 v^y 1 1 \xy 2 1 v^y 1 \x 

- (i) S S 2 * (;) " ^ (;) B ! 2> (I 
_ B «(i) Bf (I) Bf (i)- 2B f(i) B p(i 

- 2B f(I) B p.(i)- M f(I)rf(i 
^^P(i) + | B P>(i) 4 + iw 2 -l B «(i) B f(i 



(i) 

2 



-<>(I) s ™(I) +2BS ..(I) s f(I) + ^(I) + i s f(I 

- 2B f(i) + 4 B r(i) +2S f(I) + 4 B <«>(i) + 6 B f(i 

+3B f(i)-e B <"»(i)- 3B r(i) + ^'G)-i B f»(i' 
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+ ^ l5) G)^ l6) G)-^ l8) G)-^ 3) G) 31 - 2 

+ l^B? [I) log 2 + \^Bf ] 0) log 2 + Uv\oi2o(x) + 2,vr^ 6) ) a(x) 
+ ^ 3 log2a(x) - \i^a{x) - \b? Q) C 3 - §B< S > Q) Cs 



6Li4 



1\ 37^ 4 



2/ 1440 

G. Expression of HPL's in terms of the spanning set 

In this appendix we present the results for expressing all HPL's up to weight four in terms 
of the spanning set {B^}. We restrict ourselves to giving the expression for a minimal set 
of HPL's out of which all other cases can be obtained via shuffle relations 9 . 

G.l Results for weight two 

H(-l, 1; x) = log 2 log(l - x) - log(l - x) log(l + x) - \ log 2 2 - Li 2 (^^) + " 

#(0,-1; x) = -Li 2 (-x), 
H(0,l;x) = Li 2 (x). 

G.2 Results for weight three 

1 — x\ , ^ N 3, 2 , 



2 J 12 ' 



H(-l, 1, -1; x) = -Li 2 I — J log(l + x) - - log^ 2 log(l + x) 

+ log21og 2 (l +x) - log(l - x)log 2 (l + x) + log21og(l - x)log(l + x) 



+ \^ log(l + x) + I log 3 2 - ^ 2 log 2 - 2Li 3 (1±^) + ^ , 



1, 1; x) = Li 2 (j-^j Iog(l - x) - ^ log2 log 2 (l - x) + ^ log 2 (l - x) log(l + x) 

1, 3n 1 2 , „ fl — x\ 7(3 

+ - log 3 2 vr 2 log 2 - Li 3 + — , 

6 5 12 & 6 \ 2 J 8 ' 

11 7T 2 

#(0,-1, -l;x) = -Li 2 (-x)log(l + x) + -log 3 (l + x) - - log x log 2 (l + x) - — log(l + x) 

- Li3 (iT^) +C3 ' 

#(0, -1, 1; x) = Li 2 (-x) log(l -x)-^ log 3 (l - x) - ^ log 2 21og(l - x) + ^ log21og 2 (l - x) 

+ X - log 2 (l - x) logx - ^ log(l - x) + X - log 3 2 - ^ log 2 - Li 3 

f 2x \ 1 

+ Li 3 (l - x) - Li 3 (-x) + Li 3 (x) + Li 3 ( — - y ) - gCs , 



9 A set of text files containing the expressions for all HPL's up to weight four (for x € [0, 1]) in Mathe- 
matica is included in the arXiv distribution. 
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#(0,0,-1; a;) = -Li 3 (-x) , 
tf(0,0,l;x) = Li 3 (x), 

H(0, 1, -1; x) = Li 2 (x) log(l + x) + - log 3 (l - x) - - log 2 2 log(l + x) - - log 2 log 2 (l - x) 

— ^ log(l — x) log 2 (l + x) — ^ log 2 (l — x) log x + log 2 log(l — x) log(l + x) 

7T 2 7T 2 1 

+ — log(l - x) + log(l - x) log x log(l + x) + — log(l + x) + - log 3 2 

it 2 ( 2x \ ( 1 

- - log 2 + Li 3 (-x) - Li 3 (x) - Li 3 [— x j + Li 3 

- Li3 (rrf)- Li3 ( i ^) + ^' 

1 7T 2 

H(0, 1, 1; x) = -Li 2 (x) log(l - x) - - logxlog 2 (l - x) + — log(l - x) - Li 3 (l - x) + Cs 
G.3 Results for weight four 

H(-l, 1, -1, -1; x) = -hi 2 log 2 (l + x) - 2Li 3 (^±^j log(l + x) 

- ^Cs log(l + x )~\ lo S 3 2 M 1 + x) + ^ log 2 log 3 (l + x) 

- - log(l - x) log 3 (l +x)-- log 2 2 log 2 (l + x) 

+ X - log 2 log(l - x) log 2 (l + x) + ^— log 2 (l + x) + ^ log 2 log(l + x) 
+ 3Li 4 (l^)-3Li 4 (i 

tf(-l,l,-l,l;x) = ^Li 2 log 2 2-Li 2 ^I_^log21og(l-x) 

+ Li 2 log(l - x) log(l + x) + 2Li 3 log(l - x) 

- 2C 3 log(l - x) + ^Cs log(l + x) + -L log 4 (l + x) - I log 3 2 log(l - x) 

+ \ log 3 2 log(l + x) - \ log(l - x) log 3 (l + x) + log 2 2 log 2 (l - x) 

1 3 

- - log 2 2 log 2 (l + x) + - log 2 2 log(l - x) log(l + x) 

7T 2 

- 21og21og 2 (l - x)log(l + x) + log 2 (l - x)log 2 (l + x) + — log 2 (l + x) 

6 

5 7T 2 5 

+ ^g 2 log(l - X) - — log 2 log(l + X) - — 7T 2 log(l - x) log(l + x) 



1 



+ ; log 4 2 



^ 2 2 1 fl-xV 7T 2 T /1-X\ 



,X-1\ A+x\ 11 4 

+ 2L 4 - 2L 4 + vr 4 , 

1 x + 1 / V 2 J 480 ' 
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H(Q, -1, -1, -1; x) = ~U 2 (-x) log 2 (l + x) - Li 3 (A^) lo g(l + x ) + \ lo g 4 (! + *) 

1 7T 2 / 1 \ 7T 4 

- - logxlog 3 (l + x) - — log 2 (l + X)~ Li 4 ( — J + — , 

(1 2x \ 1 / 1 — x \ 
2' xTl) + 2 Li2 ( ~2~ ) l0g2(1 + X) 

- -Li 2 (-x) log 2 (l + x) + ^Li 2 (x) log 2 (l + x) + -Li 2 (-x) log 2 2 
+ Li 2 (-x) log(l - x) log(l + x) - 2Li 3 (x) log(l + x) 

- 2Lis (Si) iog(i + x) - 2Lis (tt£) iog(i + x) 

- Li 2 (-x) log 2 log(l - x) + Li 3 f lo g(! - x) - ^Cs log(l + x) 

- ^Cslog(l-x) + jhog 4 (l + x) - ^log21og 3 (l + x) 

2 1 

log(l - x) log 3 (l + x) - logxlog 3 (l + x) log 3 21og(l + x) 

3 2 

1 3 
+ - log 3 (l - x) log(l + x) + log 2 2 log 2 (l + x) + - log 2 log(l - x) log 2 (l + x) 

z 

5 

+ 21og(l - x) log x log 2 (1 + x) - — vr 2 log 2 (l + x) 

- log21og 2 (l — x) log(l + x) — log 2 (l — x) log x log (1 + x) 

+ y log 2 log(l + x) + ^vr 2 log(l - x) log(l + x) - ^Li 4 (l - x 2 ) 

1 /l — x\ 7T 2 

- -Li 2 (-x) 2 + Li 2 ( — ^— J Li 2 (-x) - — Li 2 (-x) + Li 4 (l - x) 
+ Li 4 (-x) + Li 4 (x) + hi, (^p) " 2 Li * (^r) + 2 Li4 

+ 2 Li4 (^)-2 Li4 (^) + 3 Li4 (l±^)-3 Li4 (i) + ^, 

(1 2x \ 1 / 1 — x\ 
2' xTlJ " 2 Li2 V~J l0g2(1 + X) 

+ *Li 2 (-x) log 2 (l + x) - *Li 2 (x) log 2 (l + x) - *Li 2 (-x) log 2 2 

- Li 3 [^Y^j lo g(! + x ) + Li 3(! " z) lo g(! + a:) - Lis(-aO log(l + x) 
+ 3Li 3 (x) log(l + x) + 3Li 3 (Si) lo g(! + x ) 

+ 2Li 3 ) lo g( X + x ) + Li 2(-^) log21og(l - x) + yCs log(l + x) 

- ^ log 4 (l + x) + 2 log 2 log 3 (l + x) + i log(l - x) log 3 (l + x) 

7 7 1 

+ -logxlog 3 (l + x) + - log 3 2 log(l + x) - -log 3 (l - x)log(l + x) 

7 3 

- -log 2 21og 2 (l + x) - -log21og(l -x)log 2 (l + x) 
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17 

- 21og(l - x)logxlog 2 (l + x) + — vr 2 log 2 (l + x) 

3 1 
+ - log 2 log 2 (l - x) log(l + x) - - log 2 2 log(l - x) log(l + x) 

3 7 
+ - log 2 (l — x) logxlog(l + x) — — vr 2 log 2 log (1 + x) 

- ^vr 2 log(l - x) log(l + x) + ^Li 2 (-x) 2 - Li 2 f 1 ^) Li 2 (-x) 

71-2 T- / X 1 T- / X 3 T- / X 3 T- ^ ^ 

+ tt: Li 2(-^) - o Ll 4(-a;) - o Ll 4(a;) - T Ll 4 



12 " v y 2 " v ; 2 * w 4 + 



^2 



- ^ (lix) - ^ (x^l) + ^ (xfl) - 6Li4 A + X 

+6Li4 G) + ^ 



(0, -1, 1, 1; x) = ~Li 2 (-x) log 2 (l - x) + Li 3 (^) lo g(l " " Li s(l " a) log(l - x) 

+ Li 3 (-x) log(l - x) - Li 3 (x) log(l - x) - Li 3 ^;^y^ log(l - x) 

7 19 23 1 

+ -( 3 log(l - x) + — log 4 (l - x) + — log 4 (l + x) - - log 2 log 3 (l - x) 

7 1 1 

- ^logxlog 3 (l - x) - — log(l + x)log 3 (l - x) - — log 3 (l + x)log(l - x) 

- Jlog21og 3 (l + x) - ^logxlog 3 (l + x) - Jlog 3 21og(l + x) 

o 4i o 

+ - log 2 21og 2 (l -x)-— log 2 (l + x) log 2 (l — x) + - logxlog(l + x) log 2 (l — x) 

3 1 1 

+ — vr 2 log 2 (l - x) + -logxlog 2 (l + x) log(l - x) + -log 2 21og 2 (l + x) 

13 7T 2 7T 2 

~ 48 71-2 1 ° g2( - 1 + x ) ~ 24 log ^ + ^ log ^ ~ x ) + y 1 °g 21 °g( 1 + x ) 
+ | Li 4 (1 - * 2 ) - ^Li 4 (^y) - Li 4 - Li 4 (l - x) 

1 T , x 1 T / x / x \ ( 2x \ 1 / 4x 

~ 2M-X) + ^(x) + 2L 14 ^— j - L l4 (— [ J + Z L 14 ^ 



(x + iy 



+ 2Li , ( — — ) + 2Li 4 I — — ) - 2Li 4 I ) - 2Li , 

1 + x / V x + 1 / V x + 1 



(^) 



— Li 4 



! \ 7T 4 

" 72 ' 



H(0, 0, -1, -1; x) = -Li 3 (-x) log(l + x) + ( 3 log(l + x) + — log 4 (l + x) - - logx log 3 (l + x) 

vr 2 , 21- - - • ' s ^ . / 1 \ T . / X \ 7T 4 



- - log 2 (l + x) + Li 4 (-x) + Li 4 [Y^ x )+ u ^y x + 1 j 90 , 

3 123 
ff(0, 0, -1, 1; x) = Li 3 (-x) log(l - x) + -( 3 log(l - x) + — log 4 (l - x) + - log 4 (l + x) 

- ^logxlog 3 (l - x) - ^log(l + x)log 3 (l - x) - ^-log 3 (l + x)log(l - x) 
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- i log21og 3 (l + x) - ^-logxlog 3 (l + x) - Jlog 3 21og(l + x) 

'-i o 
11 7T 2 

- — log 2 (l + x)log 2 (l -x) + -logxlog(l + x)log 2 (l - x) + — log 2 (l - x) 
16 4 Id 

1 113 

+ -logxlog 2 (l + x)log(l - x) + -log 2 21og 2 (l + x) - —ir 2 log 2 (l + x) 
4 1 4o 

- ^ log(l + x) log(l - x) + y log 2 log(l + x) + ^Li 4 (1 - x 2 ) 

I f x 2 \ 3 1 

- |Li4 ( x2 _ 1 J - Li 4 (l - x) - -Li 4 (-x) + -Li 4 (x) + Li 4 

+ 7Li 4 ( -, ^tto ) + 2Li 4 f — — ) + 2Li 4 f — ^— ] - 2Li 4 ' 



4 V( x + 1 ) 2 / Vi + ^y V^ + v v^ + i 

tf(0,0,0,-l;x) = -Li 4 (-x) , 
#(0,0,0, l;x) =Li 4 (x), 

(0, 0, 1, -1; x) = Li 3 (x) log(l + x) + -( 3 log(l + x) - - log 4 (l + x) + - log 2 log 3 (l + x) 

+ ^logxlog 3 (l + x) + ^log 3 21og(l +x) - ^log 2 21og 2 (l +x) + ^rlog 2 (l + x) 
fa 3 I 6 

7T 2 1 3 1 / 4x \ / 1 

- — log 2 log(l + x) + -Li 4 (-x) - -Li 4 (x) - -Li 4 — ^ - Li 4 



2 ^ ' 2 * K > 4 * V^ + i) 2 / U + z 

H{0, 1, -1, -1; x) = *Li 2 (x) log 2 (l + x) + Li 3 (-x) log(l + x) - Li 3 (x) log(l + x) 

- L " {Si) log<1 + x > + "» (rb) log(1 + 11 - L ' 3 (rrO log(1 + *> 

- Li 3 (^Y^j lo g(! + x ) ~ |Cs log(l + x) + 1 log 4 (l + x) - ^ log21og 3 (l + x) 

- ^log(l -x)log 3 (l + x) - ^logxlog 3 (l + x) - ^log 3 21og(l + x) 

+ * log 3 (l - x)log(l + x) + * log 2 21og 2 (l + x) +log21og(l - x)log 2 (l + x) 

+ log(l — x) logxlog 2 (l + x) — — log 2 (l + x) — - log21og 2 (l — x) log(l + x) 

8 2 

1 7T 2 7T 2 

- - log 2 (l - x) log X log (1 + x) + — log 2 log (1 + x) + — log(l - x) log(l + x) 

z fa fa 

L , , 1 T . , , 1 T . / 4x \ T . / 1 \ , ( 2x 
-Li 4 (— xj + -Li 4 (x) + -Li 4 — ^ + Li 4 - Li 



2 4V ' 2 * v y 4 ^v^ + i) 2 / Vi + ^y v^ + i 

l\ _ 7T 4 

2 J 90 ' 



+ 3Li 4 (^-3Li 4 f^- 7r4 
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47 A O 11 O 

H{0, 1, -1, l;x) = -— log 4 (l - x) + log 2 log 3 (1 -x) + — log x log 3 (1 - x) 

+ ^_ io g (i + x ) log 3 (l - x) - 1 log 2 2 log 2 (l - x) + ^ log 2 (l + x) log 2 (l - x) 

5 

- log 2 log (1 + x) log 2 (l — x) — - logxlog(l + x) log 2 (l — x) 

- (^) log 2 (l -x) + \u 2 {-x) log 2 (l - x) - l -U 2 {x) log 2 (l - x) 

- ^vr 2 log 2 (l - x) - ^ log 3 2 log(l - x) + ^ log 3 (l + x) log(l - x) 

1 7T 2 

- - logxlog 2 (l + x) log(l - x) + — log21og(l - x) 

1 7T 2 

+ log 2 Li 2 (x) log(l - x) - log(l + x)Li 2 (x) log(l - x) - Li 3 (-x) log(l - x) 
+ Li 3 (x) log(l - x) + 3Li 3 (^7y) lo g(! - x ) ~ Li 3 (i^) ~ x ) 

+ Li3 (jTlc) log(1 ~ x) + Li3 (^) log(1 ~ x) ~ f Cs log(1 " x) 

55 7 3 

- — log 4 (l + x) + log 2 log 3 (1 + x) + — log x log 3 (1 + x)-- log 2 2 log 2 (1 + x) 

29 2l 2/-, \ Li 2 (x) 2 . 1 , / 1 2x \ 

+ -vr 2 log 2 (l + x) y- - Li 2 , 2 (-1, x) + Li 2 ,2 — J 

7T 2 , „, /a , 1, o_ , , T . /1-X 



+ log 3 2 log(l + x) - — log 2 log(l + x) - - log 2 2 Li 2 (x) - Li 2 ( ) Li 2 (x) 



IT 2 



+ Li 2 (-x)Li 2 (x) + y^I^x) + 2 log(l + x)Li 3 (x) + 3Li 4 (l - x) 

1 3 T ^ x \ T /2x 
— -Li4(— x) — -Li4(x) — 2L14 + 3L14 



- 6Li 4 (1±^) - ^Li 4 (1 - x 2 ) + ^Li 4 (^) + I log(l + X) C 3 



2 J 72' 
#(0,l,0,-l;x) = -Li 2j2 (-l,x), 

H(0, 1, 0, 1; x) = 2Li 3 (x) log(l - x) - 2C 3 log(l - x) - ^ log 4 (l - x) + ^ logx log 3 (l - x) 
- ^ log 2 (l - x) + Li2 M 2 + 2Li 4 (l - x) - 2Li 4 (x) - 2Li 4 - ^ , 
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(1 2x \ 1 /l — x\ 
2'x~lJ + 2 Li2 \~) l0g2(1 



x) 



- ^Li 2 (-x) log 2 (l - x) + ^Li 2 (x) log 2 (l - x) + ^Li 2 (x) log 2 2 

- Li 2 (x) log 2 log(l - x) - 2Li 3 (^jj lo g(l - - Li s(l - x) log(l + x) 

- 2Li 3 (x) log(l + x) + ^Cs log(l - x) - ^Cs log(l + x) + ^ log 4 (l - x) 

+ | log 4 (l + x) - ? log 2 log 3 (l - x) - ^ log x log 3 (l - x) - ^ log 3 2 log(l - x) 

- - log21og 3 (l + x) - - logxlog 3 (l + x) - - log 3 21og(l + x) 

o o 

1 7T 2 3 

+ - log 2 2 log 2 (1 - x) + — log 2 (l - x) + log 2 2 log 2 (1 + x) - -7T 2 log 2 (l + x) 

2 o o 

+ log 2 log(l - x) + y log 2 log(l + x) + ^Li 4 (1 - x 2 ) + ^Li 2 (x) 2 

(1 — X \ 7T^ / 1 — X 

— ^— j Li 2 (x) - Li 2 (-x)Li 2 (x) - — Li 2 (x) + Li 4 ( — — 

( 2x \ 1 / 4x 

- 2Li 4 (l - x) + Li 4 (-x) + Li 4 (x) - 2Li 4 + -Li 4 



x-lj 2 * + 
, ( , ' ^ - f| ^ + I " ]) + 2Li 4 / ■'■ 



i + xy 2 \i + xj 2 V^ + i/ V^ + i 

H(0, 1, 1, 1; x) = -Li 2 (x) log 2 (l - x) + Li 3 (l - x) log(l - x) + - logx log 3 (l - x) 

2 4 

-^log 2 (l-x)-Li 4 (l-x) + ^, 
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